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I.  Introduction 

The  probabilistic  concepts  were  first  introduced  in  the 
1940's  (e.  g.  Ref.  1)  in  dealing  with  the  assessment  and  assur- 
rance  of  structural  safety,  and  have  been  developed  over  the 
last  three  decades  into  what  is  currently  known  as  the  struc- 
tural reliability  analysis  (e.  g.  Ref.  2).  Recently,  attempts 
are  being  made  in  various  disciplines  of  engineering  to  imple- 
ment these  concepts  in  design  situations  in  the  form  of 
probabilistic  structural  design  by  devising  reliability-based 
design  criteria  (e.  g.  Refs.  3 and  4). 

The  emphasis  of  the  structural  reliability  analysis  has 
been  placed  in  the  estimation  of  the  structural  safety  in 
terms  of  the  reliability  which  is  defined  as  the  probability 
that  a structure  subjected  to  loads  and/or  other  adverse  envi- 
ronments will  perform  its  specified  mission  v/ithout  failure. 

In  the  classical  approach,  this  probability  is  estimated  under 
the  assumption  that  all  possible  failure  mechanisms  under  the 
projected  operational  conditions  and  all  pertinent  parameters 
having  significant  effects  on  these  failure  mechanisms  are 
known  and  at  the  same  time,  probabilistic  characteristics 
(such  as  probability  distributions)  of  all  these  parameters 
are  also  known.  In  some  cases,  the  failure  mechanism  itself 
might  exhibit  intrinsic  random  characteristics  (e.  g.,  fatigue 
or  crack  initiation  process) . Every  effort  should  be  made  to 
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improve  the  confidence  of  such  an  estimate  making  use  of  all 
available  probabilistic  and  statistical  techniques.  The 
level  of  sophistication  of  these  techniques  should  not  only 
reflect  the  quality  and  quantity  of  the  information  available 
but  also  be  consistent  with  the  present  state  of  the  art  of 
(a)  structural  or  stress  analysis,  (b)  failure  analysis,  parti- 
cularly techniques  of  predicting  various  modes  of  structural 
failure,  (c)  environmental  and  load  analysis  including  pro- 
jection techniques  of  future  operational  conditions  and  (d) 
inspection  techniques  and  repair  capabilities.  This  clearly 
indicates  that  the  confidence  of  the  reliability  estimate  is 
simply  a reflection  of  the  accuracy  of  the  information  and 
the  state  of  the  art  of  the  engineering  science  pertinent  to 
the  reliability  and  maintainability  problem  at  hand. 

The  current  general  consensus  seems  to  be  such  that  the 
reliability  analysis  and  the  probabilistic  design  based  on  it 
are  referred  to  as  classical  unless  the  following  conditions 
are  met:  (1)  An  active  recognition  is  given  to  the  problem  of 

the  statistical  confidence  to  the  extent  that  an  effort  is  made 
to  estimate  such  a confidence  with  possible  use  of  the  Bayesian 
concept  and  (2)  the  problem  of  structural  safety  is  considered 
in  a wider  perspective  in  relation  to  other  significant  factors, 
such  as  economic  considerations.  Hence,  the  reliability-based 
design  criteria  in  modern  context  should  be  developed  from  the 
view  point  of  systems  analysis  with  an  emphasis  on  a systematic 
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and  consistent  use  of  analytical  tools  as  well  as  empirical 
data  available  at  the  present  time  for  the  ultimate  purpose  of 
making  design  and  management  decisions. 

In  spite  of  the  difficulties  often  associated  with  the 
structural  reliability  analysis  such  as  due  to  system  complexity, 
insufficiency  of  data,  limitation  in  the  current  technical 
knowledge,  etc.,  the  techniques  typical  of  the  reliability 
analysis  play  a central  role  in  the  development  of  the  relia- 
bility-based design  criteria:  With  the  aid  of  these  techniques, 

a definitely  better  understanding  will  be  achieved  on  (a)  the 
effect  of  increasing  the  value  of  significant  design  and 
management  parameters  (such  as  safety  and  load  factors,  stress 
allowables  and  inspection  period)  on  the  structural  reliability, 
(b)  the  sensitivity  of  the  structural  reliability  to  other 
design  and  management  parameters,  (c)  the  additional  informa- 
tion needed  to  improve  accuracy  of  the  reliability  estimate 
and  (d)  from  this  understanding,  a consistent  and  systematic 
way  of  making  design  and  management  decisions  will  emerge. 

Items  (a)  and  (b)  will  indicate  the  impact  of  design  and  mana- 
gement decisions  on  the  reliability  level  while  Item(c)  the 
cost  of  improving  the  confidence  in  the  reliability  estimate. 
Also,  sensible  and  well  disciplined  use  of  subjective  engineering 
judgement  in  the  Bayesian  framework  is  considered  potentially 
beneficial  in  certain  decision  situations  to  bring  about  the 
compromise  required  and  even  desirable  for  reasonable  blending 


3 


between  the  analytical  rigor  and  the  availability  of  pertinent 
information . 

The  purpose  to  be  achieved  by  the  design  of  aircraft 
structures  in  accordance  with  specified  reliability-based 
criteria  is  to  ensure  a quantitatively  specifiable,  optimal 
level  of  reliability  for  a fleet  of  estimated  order  of  magni- 
tude of  airplanes  intended  to  be  operated  within  the  frame  of 
an  adequately  defined  spectrum  of  missions  in  the  course  of  the 
anticipated  service  life.  The  level  of  reliability  must  be 
uniquely  related  to  the  expected  performance  of  the  aircraft 
structure  in  terms  of  the  inspection  and  repair  policy  as  well 
as  in  terms  of  the  variability  of  the  missions,  of  the  material 
performance  and  of  the  environment.  Any  significant  change  in 
these  features  should  clearly  be  reflected  in  the  change  of  the 
associated  reliability  level  so  that,  as  the  need  arises,  fleet 
management  decisions  can  be  based  on  the  rational  assessment  of 
the  consequences  of  such  changes.  Thus,  it  is  clear  that  a 
computational  routine  should  be  developed  which  would  permit 
the  quantitative  assessment  of  the  effect  on  the  reliability 
level  at  any  time  during  the  service  life,  of  any  significant 
change  in  the  conditions  under  which  the  principal  design  and 
management  decisions  have  been  derived.  In  fact,  the  establish- 
ment of  such  a routine  should  be  a part  of  the  major  effort 
in  developing  the  reliability-based  design  criteria.  Such  a 
routine  is  also  used  to  identify  those  design  and  management 
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parameters  whose  variations  produce  no  significant  changes  on 
the  reliability  level. 

In  order  to  achieve  the  purpose  outlined  above,  the 
procedures  of  probability-based  structural  design  (for  example, 
as  summarized  in  Ref.  5)  must  be  developed  considering  crack 
initiation  and  propagation  processes,  residual  strength,  ins- 
pection policy,  crack  detection  capability,  proof  test,  etc. 
so  as  to  combine  all  significant  aspects  of  aircraft  design 
that  reflect  not  only  the  selected  specific  structural  scheme 
with  its  associated  potential  failure  mechanism  but  also  the 
anticipated  inspection  requirements. 

To  achieve  the  requirement  for  actual  demonstration  and 
compliance  that  the  specified  reliability-level  has  been 
attained,  order-statistical  procedures  as  well  as  fracture- 
mechanics  considerations  must  be  combined  with  the  procedures 
referred  to  above  into  an  effective  operational  demonstration 
of  the  fleet  reliability  level,  for  example,  on  the  basis  of 
the  observation  of  the  time  to  the  initiation  of  a crack  of 
critical  or  propagating  size  in  the  first  and  possibly  second 
unit  of  the  considered  fleet. 

The  general  framework  for  the  development  of  reliability- 
based  criteria  for  aircraft,  however,  will  necessarily  emerge 
as  a compromise  between  the  application  of  such  basic  theoreti- 
cal procedures  of  probabilistic  structural  design  and  the 
requirements  of  ready  implementation  of  such  criteria.  These 
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requirements  should  provide  a format  that  makes  it  reasonably 
simple  to  translate  design-procedures  and  design-decision 
processes  currently  used  in  the  Air  Force  and  in  the  airframe 
industry  into  formally  not  too  dissimilar  procedures  and  pro- 
cesses which  reflect,  however,  the  new  probabilistic  concept 
of  the  engineering  reality.  Thus,  in  order  to  assess  the 
reliability  of  an  existing  fleet  of  airplanes  as  well  as  to 
set  the  criteria  for  the  design  of  an  airframe  for  a unit  of 
a future  fleet,  the  statistical  distribution  of  each  of  the 
significant  influencing  factors,  such  as  service  loading,  struc- 
tural performance  parameters  of  the  material  as  well  as  of  the 
fabricated  structure,  environmental  conditions,  inspection  and 
repair  procedures  must  be  adequately  characterized. 

In  this  regard,  the  present  circumstances  are  much  more 
favorable,  although  far  from  being  perfect  yet,  because  the 
immediate  availability  of  relevant  data  on  these  factors  has, 
in  recent  year,  systematically  increased  as  exemplified  by  Refs. 
6-15  on  the  gust  and  maneuver  loading,  Ref.  16  on  the 

ultimate  strength  of  aircraft  structures  and  Ref.  17  on  the 
fatigue  crack  initiation.  Also,  encouraging  is  the  fact  that 
the  use  of  such  information  either  for  the  purpose  of  relia- 
bility estimation  or  for  the  purpose  of  developing  reliability- 
based  design  and  management  criteria  has  recently  begun  to 
appear  in  the  literature:  for  example.  Refs.  18  - 22  on  inspec- 
tion, crack  detection  and  repair  of  fatigue  sensitive  structures, 
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Refs.  24  and  25  describing  a possible  statistical  approach  useful 
for  reliability  demonstration,  Refs.  25  and  26  on  the  damage  to- 
lerant requirements  and  Ref.  27  on  the  overall  reliability  study 
involving  crack  initiation  and  propagation  processes  under 
random  loading  conditions  as  well  as  optimum  inspection  policy. 

The  analytical  reliability  analysis  of  aircraft  structures 
developed  in  this  study  combines  the  material  selection  proce- 
dures and  the  material  allowables,  statistical  nature  of  loads, 
fatigue  crack  initiation  and  initial  crack  distribution  models, 
crack  propagation  and  strength  degradation,  full-scale  testing 
proof  load  testing,  and  various  crack  inspection  techniques. 

The  material  selection  is  made  on  the  primary  requirement 
that  the  aircraft  design  result  in  a light-weight  and  cost- 
effective  yet  durable  structure  with  low  structural  maintenance 
needs.  Recent  studies,  quite  often,  indicate  that  high  static 
strength  metals  tend  to  exhibit  relatively  inferior  fatigue 
performance  (not  only  the  average  trend  but  also  the  statistical 
scatter)  compared  with  medium  static  strength  metals.  This 
should  be  reflected  in  the  scatter  factor  in  the  durability 
analysis  (Ref.  28).  a similar  trend  appears  to  exist  in  the 
fracture  performance  of  metals  as  exemplified  by  a diagram 
showing  the  relationship  between  the  crack  length  parameter 
and  the  structural  efficiency  parameter  (Ref.  25).  The  material 
allowables  are  currently  based  on  relevant  MIL-HDBKs  as  well  as 
on  the  experiments  performed  to  augment  the  data  base  particularly 
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on  fatigue  and  fracture  characteristics.  The  allowables  are 
basically  responsible  both  for  the  size  of  major  structural 
components  and  for  the  durability  and  damage  tolerant  perfor- 
mance of  the  elements  thus  having  a direct  impact  on  the  risk 
of  failure. 

The  random  loading  on  aircraft  structures  considered 
herein  is  a flight-by-flight  loading  consisting  of  ground 
loads,  gust  loads  due  to  clear  air  turbulence  (including  ope- 
rational gust)  and  thunderstorms,  ground-air-ground  loads  and 
maneuver  loads.  These  loads  are  modeled  either  as  composite 
Gaussian  random  processes  consisting  of  a sequence  of  Gaussian 
patches,  or  as  single  Gaussian  processes.  Experimental  infor- 
mation tends  to  indicate  that  most  of  the  gust  loading  and  the 
maneuver  loading  for  cargo  or  transport  type  airplanes  can  be 
modeled  by  composite  Gaussian  processes  while  the  maneuver 
loading  for  fighter  type  airplanes  can  be  approximated  more 
closely  by  single  Gaussian  processes  (Refs.  6,  8 - 12). 

The  random  characteristics  of  crack  initiation  and  growth 
are  treated  in  two  different  ways.  In  one  approach,  the  time 
to  crack  initiation  is  considered  as  a random  variable  having 
a two-parameter  Weibull  distribution.  In  an  alternative  and 
newly  proposed  approach,  the  initial  crack  size  in  the  struc- 
ture is  considered  as  a random  variable.  The  analysis  method 
used  in  this  approach  does  not  depend  on  the  particular  form 
of  the  distribution  of  initial  crack  size.  For  useful  com- 
parison of  these  approaches,  however,  physically  compatible 
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distributions  of  time  to  crack  initiation  and  of  initial 
crack  size  must  be  provided.  For  this  reason,  the  probability 
density  function  of  the  initial  crack  size  is  determined  by 
a compatible  transformation  of  the  probability  density  func- 
tion of  the  time  to  crack  initiation.  Then,  aircraft  relia- 
bility study  is  performed  using  both  of  these  approaches. 

Both  approaches  use  fracture  mechanics  method  to  estimate 
propagated  crack  size  and  corresponding  residual  strength  under 
random  service  loading  conditions.  For  this  purpose,  the 
statistics  of  rise  and  fall  of  the  random  stress  process  need 
to  be  evaluated.  This  is  achieved  by  utilizing  the  approximate 
analytical  methods  given  in  Refs.  29  and  30  or  digital  simula- 
tion techniques  (Refs.  31  - 33).  The  residual  strength  is 
assumed  to  be  related  either  to  the  existing  crack  size  and 
fracture  toughness  by  the  Griff ith-Irwin  equation  in  case  of 
slow  crack  growth  design  while  it  is  obtained  as  a function  -of 
crack  size  and  other  parameters  of  redundancy  in  case  of 
redundant  design.  Once  the  residual  strength  is  established, 
the  failure  rate  is  evaluated  as  the  rate  of  upcrossing  of 
the  residual  strength  by  the  random  stress  process.  The 
probability  of  aircraft  failure  is  then  computed  on  the  basis 
of  this  failure  rate  taking  the  inspection  procedures  into 
consideration. 

Because  only  one  full-scale  model  can  usually  be  used 
for  static  (to  verify  the  design  ultimate),  fatigue,  and  damage 
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tolerance  test,  the  statistical  significance  of  such  tests 
should  be  carefully  examined  by  the  mathematical  statistics 
method  together  with  the  aid  of  the  test  performed  on  labora- 
tory specimens  and  structural  components.  Also,  it  is  highly 
important  to  clarify  the  relations  between  variability  of  the 
fracture  control  and  fatigue  parameters  observed  in  material 
tests  on  specimens  of  various  sizes  and  thickness,  and  the 
expected  scatter  for  full-scale  structural  parts  or  complete 
structures,  as  influenced  by  the  material,  the  fabrication  and 
assembly  and  the  structural  configuration  itself. 

In  the  present  study,  the  effect  of  proof  test  is  investi- 
gated by  appropriately  truncating  the  probability  density  of 
the  initial  crack  distribution. 

Investigation  of  the  impact  of  inspection  procedures  on 
the  aircraft  reliability  is  one  of  the  central  themes  of  this 
study.  The  purpose  of  inspection  is  obviously  to  detect 
fatigue  cracks  in  critical  structural  details  so  that  the  componei 
can  be  repaired  or  replaced  thus  increasing  aircraft  reliability. 
The  probability  of  crack  detection  depends  on  the  particular 
inspection  technique  used  and  it  is  usually  an  increasing 
function  of  the  crack  size.  Two  basic  types  of  inspection 
techniques  are  included  in  this  study:  (1)  Periodic  base  or 

depot  (rigorous)  inspections  utilizing  a particular  NDI  tech- 
nique and  (2)  cursory  inspections  to  detect  relatively  large 
crack  sizes  performed  periodically  between  the  rigorous 
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inspections.  The  important  factors  in  this  connection  are 
the  inspection  interval  and  the  magnification  factor  associated 
with  it  as  specified  in  MIL--A-S3444 . Determination  of  the  opti- 
mum inspection  interval  balancing  the  inspection  cost  and  the 
resulting  increase  in  aircraft  reliability  is  also  an  important 
item  to  be  considered. 

Numerical  results  are  presented  to  show  the  effect  of 
inspection  frequency,  type  of  inspection,  choice  of  crack 
propagation  model,  proof  test,  and  the  type  of  loading. 
Sensitivity  studies  are  also  performed  to  identify  and  examine 
those  parameters  to  which  the  aircraft  reliability  is  most 
sensitive . 
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II.  Random  Loading 


The  operational  flight-by-flight  service  loads  produce 
"component  stress  processes"  such  as  S^Ct),  SgC(t),  SgT(t),  Sm(t) 
and  2 (t)  due  respectively  to  ground  load,  clear  air  turbulence, 
thunderstorm  turbulence,  maneuver  and  ground-air-ground  load  in 
addition  to  the  special  loads  due  to  refueling,  terrain  follow- 
ing, etc.  (Figure  1).  The  component  process  Sj^t)  due  to  taxiing 
produces  mainly  stresses  opposite  in  sign  to  those  under  flight 
conditions  and  therefore  have  an  effect  on  the  fatigue  life. 
However,  the  effect  is  believed  to  be  insignificant  in  comparison 
to  that  due  to  other  loadings  and  is  disregarded  in  this  study. 

The  process  2 (t)  due  to  ground-air-ground  load  can  be  con- 
sidered as  a random  variable.  The  magnitude  of  this  loading  is 
large  and  it  exerts  a significant  effect  on  the  total  fatigue 
life  of  the  aircraft. 

Flight  data  tend  to  indicate  that  the  component  process  Sg(t) 

- either  Sg{-.(t)  or  SgT(t)  due  to  gust  can  successfully  be 
modeled  as  "composite  Gaussian  random  process"  composed  of  indivi- 
dual Gaussian  patches  with  different  standard  deviations,  which 
are  treated  as  random  variables.  A composite  process  model  results 
in  an  exceedance  of  the  type  exp[-x].  Figure  2 indicates  that 
such  composite  model  can  also  be  applied  in  approximation  to  maneu- 
ver loads  for  transports  and  possibly  for  bombers,  while  the 
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exceedance  of  a single  Gaussian  process  is  of  the  type  exp[vx  ] 

(see  Figure  3;  note  that  the  abscissa  is  in  x ) and  fits  well  to 
fighter  maneuver  loads. 


In  this  study,  options  are  provided  to  treat  the  component  pro- 
cesses Sg(t)  and  S^(t)  either  as  composite  Gaussian  processes 
or  as  single  Gaussian  processes. 

1-  The  Expected  Rate  of  Upcrossings 

The  expected  rate  of  upcrossings  by  the  k-th  composite 
random  stress  process  Sk(t)  can  be  defined  as 


00  (k)  + 

vk(Eo>  = Lfoi  (x)u  <v  x)dx 


ID 


(k) 

where  f (x)  are  the  probability  density  functions  of  the 

0 

random  variable  0 (standard  deviation)  and  v+(R  , x)  is  the 

o 

expected  upcrossing  rate  of  the  threshold  Rq (ultimate  strength)  by 
an  individual  stress  process  S^(t)  with  standard  deviation  x. 

For  a Gaussian  process  (t)  with  mean  XQ  and  standard  devi- 
ation a.,  v+(R  , a.)  can  be  written  (e.  g.,  Ref.  34)  as 

1 O 1 


v+(R  , 0.)  = H(l)exp{-(R  - X )2/2 0?) 
01  o o 1 
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where  is  the  characteristic  stress  frequency  (the  expected 

rate  of  zero  upcrossings)  given  by 


N 


(i) 


— (J  w2S|(w)dw/J  S|(w)dw) 

2 TT  -°°  — «> 


(3) 


in  which  S*(w)  is  the  spectral  density  of  the  stress  history 

si(t) . 
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2,  Gust  and  Maneuver  Loads  as  Composite  Gaussian 
Processes 

Gust  loads  Sa(.t)  are  modeled  as  composite  Gaussian  processes 
and  classified  into  SgC(t)  and  S gT(t)  due  respectively  to  clear 
air  turbulence  and  thunderstorm.  For  a transport  type  airplane, 
the  maneuver  loading  can  also  be  approximated  by  a composite 
Gaussian  process.  The  standard  deviations  of  individual  patches 

ck,  i = 1,  2,  (see  Fig.  1)  within  the  k-th  composite 

stress  process  are  taken  to  be  statistically  independent  and 
identically  distributed  random  variables  with  half  normal 
densities  (Refs.  8 and  9 ); 


f^ 

ai 

(x) 

= (2AaJi)J5exp{-x2/2aJ1} 

X 

> 0 

(4a) 

(2) 

f 

°i 

(x) 

= (2 Aa22)  2exp{-x2/2a22} 

X 

II V 

o 

(4b) 

f(3) 

ai 

(x) 

= (2/7ra23)  2exp{-x2/2c23) 

X 

II V 

o 

(4c) 

in  which  a , 

a „ 

and  o are  the  "intensities" 

of 

stress 

Cl 

C2 

C3 

resulting  from  clear  air  turbulence,  thunderstorm  and  maneuver, 
respectively.  It  is  implied  here  that  composite  Gaussian 
processes  associated  with  clear  air  turbulence,  thunderstorm 
and  maneuver  are  respectively  identified  by  k = 1 , 2 and  3. 
Using  Eqs . 2 and  4,  the  expected  upcrossing  rates  of  these 
composite  Gaussian  processes  are  evaluated  from  Eg.  1 as 
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Thunderstorm  (5b) 
Maneuver  (5c) 

(6) 

(7) 


in  which  S*(w)  and  S*(w)  are  the  spectral  densities  of  stress 
g m 

processes  due  to  gust  and  maneuver,  respectively.  It  is 

important  to  note  that  Eqs.  5a,  5b  and  6 are  derived  under  the 

assumption  that  the  analytical  forms  of  S£(w)  are  identical 

irrespective  of  index  i and  therefore  they  can  be  written  as 

S*(w)  multiplied  by  a constant  depending  on  i.  This  assump- 

g 

tion  also  implies  that  the  spectral  shapes  of  individual 
Gaussian  processes  S^(t)  are  identical  whether  due  to  clear 
air  turbulence  or  thunderstorm.  The  same  assumption  is 
used  for  Eqs.  5c  and  7. 

The  expected  upcrossing  rate  under  the  combined  loading 
is  then  determined  as 


v+(V  - = hVxp,-(Ro ' V/aci} 


+ P2NgeXf><-(R0  - Xo)/ffc2)+  P,V*pMro  - Xo)/ac3>  (8) 
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in  which  p^,  and  p^  represent  the  fraction  of  clear  air 

turbulence,  thunderstorm  and  maneuver,  respectively.  It 

should  be  noted  that  p + p + p = 1 . 

*1  2 3 

3.  Gust  and  Maneuver  Loads  as  Single  Gaussian  Processes 
When  gust  and/or  maneuver  loads  are  taken  as  single 

Gaussian  processes  S (t)  , the  probability  density  function  in 

K. 

Eq.  1 can  be  replaced  by  a delta  function 
(k) 

f (x)  = S(x  - a ) (9) 

sk 

in  which  a is  the  standard  deviation  of  the  stress  corres- 
sk 

ponding  to  a single  Gaussian  process.  Then,  from  Eqs.  1 and 
9,  the  average  upcrossings  are 

+ (k) 

v,  (R  ) = N exp{ - (R  - X ) 2 /2o2  } (10) 

k o o o sk 

( k ) 

where  N is  obtained  from  Eq . 3 by  replacing  S?  (w)  by 

S*  (w)  the  spectral  density  of  S (t)  . 

K K 

4.  Combination  of  Composite  and  Single  Gaussian  Processes 
It  should  be  noted  that  with  this  approach  a segment  of 

the  random  load  can  be  treated  as  a composite  process,  while 
the  remaining  portion  can  be  taken  as  single  Gaussian  processes. 
For  example,  if  the  gust  is  considered  as  composite  processes 
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and  the  maneuvers  as  a single  Gaussian  process,  the  expected 
rate  of  upcrossing  under  the  combined  effect  is 

V+(R  ) = p,  N exp{-(R  - X )/o  } + p N exp{-(R  - X ) /a  } 

o lg  o o cl  2g  ooC2 

+ p N exp{ - (R  - X ) 2/2o2 _}  . (11) 

3 m o o hi 

where  Nra  = 1JV  is  obtained  from  Eq.  7 with  S^(w)  replaced  by 

I 

S3 (w)  (spectral  density  of  maneuver  load  process  83(1)  = Sm(t) 

idealized  as  a single  Gaussian  process) . 
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III.  Fatigue  Crack  Propagation 


1.  Crack  Propagation  under  Constant  Amplitude  Load 

The  fracture  mechanics  theory  is  applied  to  cracks  of 
detectable  size  initiated  at  certain  time  or  to  preexisting 
cracks  for  the  purpose  of  determining  their  propagating  size 
under  a stress  history.  Under  a constant  amplitude  stress 
history,  crack  propagation  can  be  represented  by  the  well 
known  power  law  formula  (e.  g.  Refs.  35  - 39) 

da  b 

— = c(AK)  AK  > AK  (12) 

dn  TH 

where  da/dn  is  the  rate  of  crack  propagation  (per  cycle) , "a" 
is  crack  size,  Ak  is  range  of  stress  intensity  factor  and  c 
and  e are  material  constants.  These  material  constants  exhibit 
statistical  variations  from  specimen  to  specimen.  However, 
within  a specimen,  it  appears,  they  can  be  treated  as  deter- 
ministic quantities.  The  values  of  these  constants  are  deter- 
mined from  the  experimental  results  such  as  those  shown  in  Figs. 
4 and  5.  The  value  of  b equal  to  4.0  has  been  used  by  many 
investigators.  It  is  also  pointed  out  that  below  the  threshold 
value  AKm„,  the  rate  of  crack  propagation  is  drastically  reduced 
as  observed  in  Figs.  4 and  5 and  is  so  treated  in  the  present 
analysis . 
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2.  Crack  Propagation  under  a Random  Loading 

The  available  experimental  and  empirical  information  on 
crack  propagation  under  random  loading  suggests  that  a modified 
form  of  the  power  law  given  in  Eq.  12  be  constructed  to  repre- 
sent the  fatigue  crack  growth  under  such  loading.  In  the 
present  study,  based  on  the  experimental  data  given  in  Figs. 

6 and  7,  it  is  assumed  that  the  crack  propagation  under  random 
loading  follows  the  law  corresponding  to  the  bilinear  relation- 


(13a) 

(13b) 


where  c,  c',  A and  X'  are  material  parameters  to  be  determined 

from  Figs. 6,  7 and  9.  The  expected  value  of  the  b-th  power 

E 

of  the  range  of  stress  intensity  factor,  AK  , can  be  related 

to  the  expected  value  of  the  b-th  power  of  rise  and  fall  of 

b , 

the  stress,  S , and  the  crack  length  a as  follows  (Ref s .29 , 37 , 38)  . 


~b  b/2 

S (ira/2) 


(14a) 


E l/b  ~E  !/b 

(AK  ) = (S  ) (fra  /2) 

W TH 


(14b) 


Eq.  14b  provides  the  definition  of  a . The  last  two  equations 

TH 

are  valid  for  a single  through  crack  in  an  infinite  plate  with 
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It  follows 


stress  intensity  factor  K given  by  K 

from  Eqs.  13  and  14  that 

da  — p-  A/b  A/ 2 A/2 

— = c ( S ) ( tt/2)  a 

dn 

da  -r-  A'/2  A'/2  A'/2 

~ = c*  ( SD)  (tt/2)  a 

dn 

The  present  formulation  is  based  on  the  existing  practice 
of  data  processing  where  the  logarithm  of  the  rate  of  crack 
propagation  under  random  loading  is  plotted  against  (AKb) 
with  the  same  value  of  b as  used  in  Eq.  12.  As  Fig.  9 illust- 
rates, however,  the  experimental  result  indicate  that  the 
plot  does  not  necessarily  result  in  a straight  line  and  that 
the  slope  of  such  a plot  is  generally  different  from  b mentioned 
above.  In  fact,  the  present  study  uses  the  bilinear  relation- 
ship with  two  distinct  slope  values  A and  A ' and  a threshold 

value  (AK^)^k. 

TH 

Solving  Eqs.  15a  and  15b,  fatigue  crack  length  at  time 

t can  be  obtained  from 
n 

1-  A/2  A/2  A/b  1-A/b  2/(2-A) 

a(t  ) = {a  + (1  - A/2)c(tt/2)  Q N t } 

no  on 

a > a (16) 

o = TH 

l-A'/2  A'/2  A'/b  l-A'/b  2/ ( 2-A  * ) 

a ( t ) = (a  + (1  - A'/2)c(tr/2)  Q N t } 

no  on 

a £ a (t)  < a (17) 

o - --  TH 


= S(TTa/2)1/2. 


; e> . i/b 
) > 

(AKb)^b 

TH 

(15a) 

,b.  1/b 

■ ) < 

(AKb) 1/b 
TH 

(15b) 
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2/  (2-A ) 


l-A/2  A/2  A/b  1-A/b 

a (t  ) = {a  +(1-  A/2)  c (tt/2)  Q N (t  - t )} 

n TH  on  TH 

a < a < a (t)  (18) 

o — TH  = 


where  a is  the  initial  crack  size,  t = t - t , t is  time 
o n o o 

to  crack  initiation  (this  is  equal  to  zero  if  preexisting 

‘ 

crack  is  considered) , Q = N S with 

o 


P = (PlN  S*  + p N Sb  + p N S + N S^)/N 
1 g gc  2 g gT  3mm  z z 0 


(19) 


N = p N + 

P„N 

+ p.. 

N + N 

(20) 

o 1 g 

2 g 

3 

m z 

2 1-A 

1 /2 

1-A 

’/2 

A’/b  1-A '/b 

A ' /2 

bTH  2_\ , {aTH 

— 

a 

o 

}/{C 

'Q  NQ  (tt/2) 

} 

(21) 

In  Eqs.  19  - 21, 

b 

“b 

b b 

s , 

gc' 

SgT' 

s , s 
m z 

are  the  expected 

values 

of  the  b-th  power  of  rise  and  fall  of  stress  histories  associ- 
ated with  clear  air  turbulence,  thunderstorm,  maneuver  and 
ground-air-ground  load,  respectively. 
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IV.  Statistics  on  Rise  and  Fall  of  a Random  Process 

As  observed  in  Sec.  Ill,  the  rate  of  crack  propagation 
is  directly  proportional  to  the  expected  value  of  the  b~th 
power  of  the  rise  and  fall  of  the  random  stress  process. 

If  the  stress  process  S (t)  associated  with  a particular 
loading  (e.  g.  clear  air  turbulence,  thunderstorm,  etc.)  is 
a single  Gaussian  with  mean  zero  and  standard  deviation  a 

Sk 

the  expected  value  can  be  written  in  the  form 

where  A is  determined  either  by  using  the  approximate 
analytical  techniques  proposed  in  Refs.  29  and  30  or  by  the 
method  of  random  process  simulation  (Refs.  31  - 33).  If  the 
stress  process  S (t)  is  a composite  process  consisting  of 

K 

patches  of  single  Gaussian  processes  S^(t)  with  mean  zero  and 

standard  deviation  o , the  expected  value  of  the  b-th  power 

i 

of  the  rise  and  fall  of  S^(t)  is  given,  similarly  to  Eq.  22, 
by 

Sb(o  ) = A ab  (23) 

i i 

and  the  expected  value  of  the  b-th  power  of  the  rise  and  fall 
of  the  composite  process  S (t)  is  given  by 

K. 


22 


b , b (k) 

S.  = / S (x)f  (x)dx 
k J o a 


(24) 


(k) 

where  f is  the  density  function  of  a.  as  introduced  in 
o 1 

Sec.  II. 

From  Eqs.  4a,  4b,  4c,  22,  23  and  24,  for  clear  air 

K 

turbulence  (k  = 1) , thunderstorm  (k  = 2)  and  maneuver  (k  = 3) 
are  obtained  as 


b 2n 

S = A { (2n  - 1) I ! }a 
gC  lb  Cl 


(b  = 2n) 


, h n 2n+l 

= A].b(2/u)  n!2  aci 


(b  = 2n  + 1)  (25) 


sgT  A2b{ (2n  " 1} ! !}ac2 


(b  = 2n) 


% n 2n+l 
= A2b(2/u)  n!2  aC2 


(b  = 2n  + 1)  (26) 


S£  - A3b((2n  - 1)U>aC3 


(b  = 2n) 


^2  n 2 n + 1 

= A (2/ir)  n!2  ac3 


(b  = 2n  + 1)  (27) 


with  (2n  -1)!!  = 1*3*5' 


(2n  - 1) 


In  Eqs.  25  - 27,  A , A^  and  A0l  represent  the  coefficient 

lb  2b  ^>b 

A in  Eq.  23  associated  with  the  stress  due  to  clear  air 
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turbulence,  thunderstorm  turbulence  and  maneuver,  respectively. 

For  the  purpose  of  evaluating  these  coefficients  by 
simulation  method,  sample  time  histories  of  S(t)  are  generated 


by  means  of  the  simulation  technique  given  in  Refs.  31  - 33 
assuming  the  normalized  spectral  density  function  shown  in 
Fig.  10  which  can  be  expressed  as 


S(w)  = — 
(1 


B ) w 
c c 


= 0 


for  B w < w < w 
c c - = c 


otherwise 


(28) 


where  6 = measure  of  bandwidth  of  S- (t)  and  m = upper  cutoff 

c i c 

frequency.  Portions  of  these  simulated  time  histories  are 
shown  in  Figs.  11-16  for  various  values  of  parameter  B . 
Then,  the  average  of  the  b-th  power  of  rise  and  fall  of  S^(t) 
is  computed  using  the  two  counting  methods  indicated  in  Fig. 
17.  In  the  first  method  (case  I) , every  one  of  rise  and  fall 
is  computed,  raised  to  the  power  b and  averaged,  while  in  the 
second  method  (Case  II)  the  rise  and  fall  are  constructed  by 
considering  only  those  alternating  peaks  and  troughs  which 
are  the  maxima  or  the  minima  in  the  time  intervals  between 
successive  zero  crossings.  The  results  are  summarized  in 
Figs.  18  - 21  for  various  values  of  b.  The  irregularity 
factors  evaluated  from  the  sample  time  histories  compare 
extremely  well  with  the  theoretical  result  as  shown  in  Table 
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1.  Also,  a comparison  between  the  values  of  A obtained  form 
simulation  and  the  ones  by  the  approximate  theoretical 
methods  (Refs.  29  and  30)  is  presented  in  Table  2.  A good 
agreement  is  observed  for  those  values  of  b considered  in 
the  theoretical  methods  under  the  conditions  of  Ref.  30 
(narrow  Band  with  = 1.0)  and  of  Ref.  29  (wide  band  with 
Bc  = 0) . The  results  for  narrow  band  cases  are  particularly 
useful  with  respect  to  the  stress  response  process  to 
turbulence,  since  such  response  process  is  expected  to  be 
narrow  band  because  of  structural  filtering  effects, 

A collection  of  many  acceleration  records  during  maneuver 
indicates  that  the  (acceleration)  load  history  exhibits  the 
characteristics  of  a wide  band  random  process.  To  represent 
the  maneuver  load  in  frequency  domain,  therefore,  the 
following  form  of  spectral  density  has  been  empirically 
suggested  for  a fighter  type  airplane  (Ref.  12  ). 

SQ 

S (w)  = oj_  < a)  < 00  (29) 

m 2tt{  1 + ((d/(dB)d}  D “ 

2 

where  d = 2.6,  = 2 tt  x 0 . 031  (rad/sec)  and  SQ  = 30  (g  -sec) 

with  <dD  indicating  the  lower  cut-off  frequency.  If  the 

maneuver  loading  is  assumed  to  be  a (single)  Gaussian  random 
process,  then  the  load  exceedance  curve  (interpreted  as  in- 
dicating level  crossing)  is  symmetric  with  respect  to  negative 
and  positive  loading.  Experimental  measurements  tend  to 
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indicate  (as  shown  in  Fig.  2)  that  such  an  assumption  may  be 
valid  for  an  airplane  of  transport  type.  However,  as  observed 
experimentally  and  also  as  specified  in  Military  standards 
(see  Fig.  3),  the  positive  and  negative  parts  of  the  load 
exceedance  curve  for  a fighter  due  to  maneuvers  are  markedly 
different.  To  account  for  this  fact,  the  following  procedure 
is  proposed:  Consider  two  Gaussian  processes  Y (t)  and  Yn(t) 
which  are  associated  with  the  symmetric  exceedance  curves 
respectively  constructed  from  positive  and  negative  parts  of 
the  (specified)  asymmetric  exceedance  curve.  Let  ap  and 
denote  the  standard  deviations  of  Y (t)  and  Y^(t) • These 
standard  deviations  can  be  evaluated  from  the  positive  and 
negative  slopes  of  the  specified  asymmetric  exceedance  curve. 
Then,  digitally  generate  Y (t)  in  the  same  manner  as  above 
using  the  spectral  density  given  in  Eq . 29  (see  Fig.  22  for 
sample  time  history)  and  construct  Y(t)  so  that 


Y ( t ) = 


Yp(t) 


rYp (t) 


Y (t)  > 0 

P 


(30) 


Y (t)  < 0 


where  r = an/a  . The  process  Y(t)  thus  constructed  produces 

the  same  exceedance  curve  as  specified.  A sample  time  history 

of  Y (t)  with  a)  = 2tt  x 0.007  rad /see  and  r = 0.3  is  shown  in  Fig.  23. 
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From  the  time  history  Y(t),  the  expected  value  of  the 
b-th  power  of  rise  and  fall,  S^,  is  calculated  using  the  coun- 
ting procedures  indicated  in  Fig.  17.  It  should  be  noted  that 
the  parameters  d,  w , in  Eq.  29  can  be  adjusted  to  repre- 

sent maneuver  loading  for  a variety  of  fighters  under  different 
flight  conditions.  The  transformation  between  the  maneuver 
loading  in  terms  of  acceleration,  Y(t),  and  stress,  Sm(t) , can 
be  performed  in  a quasi-static  manner.  This  is  due  to  the  fact 
that  maneuver  loads  are  applied  at  a relatively  slow  rate  in 
comparison  to  the  structural  dynamic  sensitivity  of  a critical 
part  in  an  aircraft.  The  energy  in  the  maneuver  spectral 
density  Sm(w)  is  concentrated  at  low  frequencies  while  the  peaks 
of  the  absolute  frequency  response  function  |H(w) | appear  at 
higher  frequencies.  Thus,  the  time  history  of  the  stress  at 
a particular  critical  point  in  the  aircraft  due  to  maneuvers  can 
be  written  as  Sm(t)  = CY (t)  where  C is  a proportionality  con- 
stant. For  example,  from  Ref.  13  (Fig.  24)  it  can  be  observed 

2 

that  C - 4KGF/mm  for  wing-root  bending  stress  for  a fighter. 

Because  of  this  proportionality,  the  coefficient  A in  the 

TT" 

equation  = Aab  (see  Eq.  22)  is  identical  with  the  coefficient 
A in  the  equation  Y*3  = Aa  where  Y^  is  the  expected  value  of  the 

Xr 

b-th  power  of  rise  and  fall  of  Y(t).  Hence,  the  values  of  A 
are  evaluated  from  time  histories  of  Y(t)  (Fig.  23)  and  plotted 
in  Fig.  25  as  functions  of  b (dotted  curves) . The  solid  curves 
in  Fig.  25  indicate  the  values  of  A appearing  in  the  equation 
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= Acr*3  and  obtained  from  time  histories  of  Y (t)  (Fig.  22)  . 

P P P 

The  dotted  and  solid  curves  in  Fig. 25  , therefore,  show  the 
difference  between  based  on  the  specified  asymmetric  excee- 
dance curve  and  on  the  corresponding  symmetric  exceedance  curve. 

It  is  pointed  out  that  if  Y (t)  is  a narrow  band,  evaluated 
* p m 

from  Y (t ) is  equal  to  evaluated  from  Y (t)  multiplied  by 

b Id 

(1  + r)  /2  . In  the  numerical  example  that  follows,  the  solid 
curve  for  Case  II  is  used  for  aircraft  of  transport  type  because 
they  tend  to  have  symmetric  exceedance  curves  as  previously 
mentioned  while  the  dotted  curve  for  Case  II  is  used  for  fighter 
aircraft . 
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V. 


Residual  Strength 


1.  Slow  Crack  Growth  Model 

When  the  crack  is  interpreted  as  initiating  rather  than 

pre-existing,  a crack  of  size  aQ  is  initiated  at  time  tQ  after 

the  aircraft  is  placed  into  operation.  Then,  the  initiated 

crack  grows  in  its  size  in  the  manner  as  described  in  Sec.  Ill, 

and  under  the  slow  crack  growth  design,  the  (initial)  ultimate 

material  strength  Ro  starts  to  decrease  as  crack  reaches  the 

initial  critical  crack  size  a.  at  time  t.  (=  t + t.  ) where 

t,  is  the  time  required  for  the  crack  to  increase  its  size 
ic 

from  a„  to  a . . 

o ic 

The  relationship  between  the  residual  strength  and  the 
crack  size  can  be  expressed  by  the  equation  of  Griff ith-Irwin 
type 

I<  = R0/Traic/2  (31a) 

K = R ( t ) /na(t  )/2  (31b) 

C m m 

where  is  usually  larger  than  and  = critical  stress  inten- 

sity factor,  a(t  ) = crack  size  at  time  tn  + t ; i.e.,  at  time  t 

m 1mm 

after  the  crack  reaches  a^c>  The  interrelationship  among  Rq, 

R(tm),  aQ,  a^  , a(tm),  etc. is  schematically  illustrated  in  Fig. 26 

while  the  relationships  among  t , t,  , t , t . and  t (t  will  be 

3 o 1 m ic  nn 

defined  later)  are  illustrated  in  Figure  27, 
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The  time  t^  is  evaluated  from  Eq.  16  or  Eq.  17  as 


t . 
xc 


2 


2 - 


l-A/2 
(a . 


A 1C 


l-A/2  A/b  1-A/b  A/2 

a )/  cQ  N (u/2) 

o o 


t . 

1C 


2 

+ 


2 


2 


A 

2 


(a 


l-A'/2 

TH 


l-A/2 
(a . 


a.  > a > a (32a) 

xc  o - TH 


l-A'/2  A'/b  l-A'/b  A'/2 

a )/  c'Q  N (tt/2) 

o o 


l-A/2  A/b  1-A/b  A/2 

aTH  )/  cQ  (it/ 2) 


a < a (32b) 

o TH 


Then,  the  expression  for  the  residual  strength  as  a function 

of  time  t becomes 
m 


R(t  ) 
m 


% 


l-A/2 


R a . /{a 
o ic'  ic 


A/2  A/b  1-A/b  l/(2-A) 

(1  - A/2)  c (tt/2)  Q N t } 

o m 


(33) 

It  is  important  to  note  that  all  the  equations  and  the 

definitions  of  parameters  appearing  in  this  subsection  are 

valid  even  when  the  crack  of  size  a is  assumed  to  be  pre- 

o 

existing  (rather  than  initiating  at  tQ)  so  long  as  t is  taken 
to  be  zero. 


2.  The  Fail-Safe  Design 
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A design  practice  to  prevent  the  reduction  of  the  residual 
strength  from  reaching  an  excessive  level  is  to  introduce  crack 
stoppers  in  the  structure,  thus  making  the  design  fail  safe. 

In  this  case,  Eqs.  31a  and  31b  are  no  longer  valid  to  represent 
the  strength  of  a cracked  structure.  The  strength  of  the  fail 
safe  structures  depends  on  a particular  design  and  the  residual 
strength  should  be  determined  by  individual  analysis  and  testing. 
Following  Refs.  27,  40  and  41,  the  residual  strength  after  tn 
flight  hours  is  assumed  to  be 

% 

a (t  ) - a 

R ( t ) = R {1  - (1  - V,  ( ) } (34) 

n ° a - a 

s o 

where  ag  = fail-safe  crack  size,  RQ^  = residual  strength  at  the 
fail-safe  crack  size  ag,  and  0 < £ < 1. 
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VI.  Failure  Rate 


A failure  of  the  structure  occurs  when  the  residual  strength 

R(t^)  is  exceeded  by  the  applied  (random)  stress.  Then,  the 

problem  is  essentially  that  of  a first-passage  probability  with 

a variable  one-sided  threshold  (Refs.  42  - 44).  The  expected 

upcrossing  rate  v+(R  ) of  the  stress  process  S,  (t)  with  respect 

k o k 

to  the  threshold  value  Rq  (initial  resisting  strength)  obtained 
in  Sec.  II  can  also  be  used  in  approximation  for  the  expected 

rate  with  the  variable  threshold  R(t  ).  The  expected  failure  rate 
hj,'(tn,  Rn)  / or  risk  function,  associated  with  this  first 
passage  problem  can  then  be  approximated  by 

h (t  , R ) = v+(R  )/m  (35) 

k n n k n c 

where  R is  written  for  R(t  ) and  the  clumpsize  M > 1 (Refs, 
n n c = 

43  and  44).  For  aircraft  structures,  the  difference  between 

the  initial  strength  R and  the  stress  X due  to  one  g load 

o o 

is  by  design  much  larger  than  the  standard  deviation  of  the 

stress  process, and  also  the  restriction  on  the  length  of 

service  life  as  well  as  the  implementation  of  the  inspection 

maintains  the  residual  strength  R(t  ) at  a higher  level  compared 

n 

with  Xq  in  terms  of  the  standard  deviation  of  the  stress  process. 
Therefore,  the  events  of  stress  excursion  beyond  the  threshold 
can  be  assumed  to  be  statistically  independent  and  hence 
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Mc  = 1 is  used  for  the  present  purposes. 

In  this  investigation,  the  residual  strength  is  also  con- 
sidered as  a random  variable  and  therefore,  the  failure  rate 
given  in  Eq.  35  with  M = 1 is  modified  as 

X 00  + 

h (t  ) = / f (x)dx  + / v,(x)f  (x)dx  (36) 

K n -00  K X 

O 

where  f (x)  is  the  probability  density  function  of  the  residual 
strength  R(tn).  In  this  equation,  the  value  of  unity  is  used 
for  simplicity  in  place  of  v*(x)  for  the  values  of  x less  than 

.K 

X . Since  this  is  the  range  of  x where  fD(x)  is  insignificantly 
small,  the  resulting  value  of  the  integration  is  rather  insensi- 
tive to  the  assumption  of  v£(x)  in  this  range  of  x. 

A further  modification  is  necessary  when  the  loads  on  the 
aircraft  can  be  controlled  by  pilots  or  by  other  means.  For 
example,  if  the  controlled  maneuver  results  in  a truncation  of 
the  load  spectrum,  the  failure  rate  is  obtained  from 


(37) 

where  mT^0  is  the  truncation  level  with  mT(>l)  being  the  trun- 
cation index.  To  determine  the  failure  rate  due  to  a combined 
effect  of  each  component  process  k,  one  can  use 
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(38) 


h,V  = J/kW 


For  a Gaussian  distribution  of  the  residual  strength  R(t  ) 
with  mean  value  Ynyo  an(3  standard  deviation  aR  = v0Yny0/  the 


probability  density  function  is 


f„(x)  = exp{  ■ 

“ /2itor 


(x  : Vo1*, 

2V 


(39) 


where  for  the  slow  crack  growth  model,  v is  obtained  from 

n 

Eq.  33  as 

2~*  _1_ 

1 2~  X/2  A/2  1-A/b  2-A 

Y =a  V(a.  + (1  - A/2 ) c (tt/2  ) Q N (t  - t.  ))  (40) 

n ic  ic  o n ic  J 


and,  for  the  fail-safe  design  model,  from  Eq.  34, 


Y 


n 


1 - (1  - K){ 


a ( t ) - a 


n 


a - a 
s o 


?.} 


(41) 


When  t = t in  Eq.  40  and  t = 0 in  Eq.  41,  y becomes  unity 
n ic  n n 

and  the  density  function  in  Eq.  39  represents  that  of  the 


initial  resisting  strength  R . 

o 


+ 


Using  Eq.  39  in  Eqs.  36  and  38  with  expression  of  v (x) 

K 

given  in  Eq.  5,  one  obtains  the  failure  rate  associated  with  a 
combined  stress  process  with  no  load  truncation  as 
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(42) 


3 n . 

h(t  ) = l %P.N.{  [1  - erf(-  1— )]  + exp[-(2f).  - r2)/2] 
n i=l  /2r . 1 


n±  ~ rl 

[1  + erf  (-  -)  ] } 

/Zr 


in  which  Nn  = N = N , N = N , n = (y  y /a  .)  - (X  /a  .) , 
i 2 g 3 mi  n o ci  o ci 

ri  = voVo/oci'  Vo  = °E  /vo'  1 =1’  2'  3- 

o 

When  each  of  the  component  loading  processes  is  assumed 
to  be  a single  Gaussian  with  a truncated  average  upcrossing 
rate,  the  failure  rate  associated  with  the  combined  stress 
is  obtained  from  Eqs.  10,  37,  38  and  39. 


(t  ) = l %P.N.{(1  -erf  ( — ) + (1  + r?)  zexp[-%( — + — ) 

n i=l  1 1 /2 F.  1 V2  V2 


n . 


ri  1 -i 

+ %(—  + — )2(i  + r2.)  ] [erf  ( — (1  + r2. ) 2) 

1 X 


Vx 

x o 


/2ni 


+ erf  { ~ -rr— :[  (m  - 1)  1 + m 

/2  /r.  + 1 T v 1 r.V 

k x i x 


o 


}]} 


(43) 


where  V = a /X  . 

x s o 

If  the  combined  stress  process  consists  of  single  and 
composite  Gaussian  processes,  the  failure  rate  can  be  expressed 
in  an  appropriate  combination  of  terms  appearing  in  Eqs.  42 
and  43.  The  computer  program  is  written  so  that  an  arbitrary 
combination  of  this  nature  can  be  accommodated. 
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VII.  Periodic  Inspections 


The  purpose  of  inspection  is  to  detect  the  fatigue  and 
pre-existing  cracks  in  the  structural  components  so  that,  before 
cracks  become  cr itical , they  can  be  replaced  by  uncracked  com- 
ponents to  endure  their  designed  initial  strength  at  least  at 
the  time  of  replacement.  In  this  study,  two  types  of  aircraft 
inspection  at  different  levels  of  sophistication  are  considered: 
(1)  "rigorous"  periodic  inspections  which  might  correspond  to 
the  depot  or  base  level  inspection  (Ref.  45)  utilizing  parti- 
cular non-destructive  inspection  (NDI)  techniques  and  (2) 
"cursory"  inspections  which  might  correspond  to  in-flight  evi- 
dent, ground  evident,  walk  around  and  special  visual  inspections 
(Ref.  45)  . 

The  probability  of  detecting  a fatigue  crack  at  a struc- 
tural detail  during  a rigorous  inspection  depends  on  the 
probability  of  inspecting  this  cracked  detail  and  the 
resolution  capability  of  the  particular  inspection  technique. 
Typical  NDI  techniques  presently  used  include  delta  scan,  shear 
wave  ultrasonic,  magnetic  particle,  X-ray,  and  magnetic  rubber 
inspection  (MRI)  methods. 

Define  as  the  probability  of  inspecting  a cracked 
detail  and  U„ (a)  as  the  probability  of  detecting  a crack  of 
size  a.  Then,  the  probability  of  detecting  a crack  during  a 
rigorous  inspection  can  be  obtained  in  general  from 


36 


where  it  was  assumed  that  and  U2(a)  are  independent.  The 
probability  of  not  detecting  a crack  F* (a)  is  obviously  equal 
to  1 - F (a)  . 

Since  at  this  time,  the  information  on  the  probability 
is  limited,  it  is  assumed  in-  this  study  that  U-^  = 1,  i.  e., 
every  critical  detail  will  be  inspected.  Based  on  the  experi- 
mental and  empirical  results,  however,  the  detection  probability 
U 2 (a)  may  be  constructed  in  the  following  fashion  (Refs.  21,  25 
and  27  ) . 

U2(a)  = 0 a < ax 

m 

= [ (a  - a1)/(a2  - a^  ] ag  a £ a2  (45) 

=1  a2  < a 

where  1/8  < m < 1/5  for  accurate  NDI  techniques,  m > 1/5  for 
more  crude  NDI  techniques,  a-^  = the  minimum  crack  size  below 
which  the  crack  cannot  be  detected  with  the  particular  detection 
technique  used,  a2  = the  maximum  crack  size  beyond  which  the 
crack  can  always  be  detected. 

In  Fig.  28  . a plot  of  U2(a)  is  shown  for  two  different 
values  of  parameter  m (=1/8  and  1)  with  a-^  = .02"  and  a2  = 
Detection  probabilities  experimentally  obtained  for  X-ray, 


.3"  . 


ultrasonic  and  dye  penetrant  techniques  are  also  indicated  in 
this  figure  (Ref.  25). 

A similar  crack  detection  probability  model  can  be  const- 
ructed for  the  cursory  inspections.  However,  since  no  definite 
experimental  information  is  presently  available,  the  following 
simple  form  is  assumed  for  the  probability; 


Fc(a)  = 0 
Fc(a)  = 1 


a < a| 

a > ai 


(46) 


where  a'  = the  crack  size  below  which  the  crack  cannot  be 
1 

detected  by  cursory  inspection  and  beyond  which  the  crack  is 

always  detected.  The  probability  of  not  detecting  a crack  by 

the  cursory  inspection  is  F* (a)  = 1 - F (a) . 

c c 
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VIII.  Probability  of  Failure 


1.  Probability  of  Failure  Based  on  Distribution  of  Time  to 
Crack  Initiation 

Following  Ref.  27  for  the  development  of  this  subsection, 

let  PQ  be  the  probability  of  failure  of  the  aircraft  within 

the  intended  service  life  T with  no  inspection.  Then,  assuming 

the  time  t to  crack  initiation  to  be  a random  variable  with 
o 

a probability  density  function  W(tQ),  the  probability  of  failure 
can  be  expressed  as 

T 

P = 1 - exp{-Th  - (T/B)a)  - f W(t)exp{-th  -H(T  -t) }dt  (47) 
o o J o 

0 

where  hQ  is  the  expected  failure  rate  corresponding  to  the 
threshold  R„  or  h - h(t  ) in  Eg.  38,  and 

t 

n 

H ( t ) = / h ( t) dt  (48) 

n 

0 

In  this  study  as  in  Ref.  27,  a two-parameter  Weibull  distri- 
bution of  the  time  to  crack  initiation  is  used  for  W(tQ); 

a t a “I  t a 

W(t  ) = '•■(•  °)  exp{-(— ■•)  } t 0 (49) 

° 8 B 8 ° ~ 

However,  when  the  minimum  life  T*  is  reliably  known,  the 
following  three  parameter  Weibull  density  can  advantageously 
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be  used  (Ref . 46 ) . 


* a"1 

0£  t — T*  Ot 

W(t  ) = -(— - -)  exp{ - [ (t  - T*)/0]  } t > T (50) 

° 0 3 ° 0=0 

in  which  a = shape  parameter  and  3 = scale  parameter. 

Assume  now  that  the  aircraft  is  subjected  to  a rigorous 
inspection  at  the  end  of  each  Tq  flight  hours  as  indicated  in 
Fig.  29.  The  probability  of  failure  P(j)  within  the  [o,  jT  ] 
interval  ( j — 1 inspections)  can  then  be  obtained  from 


P(j) 


j-1  T 

I / °q.  . (t)W[(i-l)T  +t]dt 
1=1  J ID  ° 


(51) 


j = 2 , 3 , ... 
i = 1,  2, , j-1 


where 


(1)  D_1 

q . . (t)  = F [a (T  -t) ] C . . ( t)  +{  n F*[a(kT  -t) ] }V.  . (t) 

ID  o x]  k=  2 o 13 


D~i  k_1  m 

+ 6.  . „ l { n F* [a (mT  -t) ] } F [a (kT  -t) ] C . ^ (t)  (52) 

k=2  m=l  ° o 


± a a 

P*  = exp{-[ (j-l)To/3l  }-  exp{-[jTQ/3]  - jT^} 

T 

•Lo 

- f W[ (j-1) T +t]exp{-h  [ ( j -1 ) T +t]-H(T  -t) }dt  (53) 

1 O O O o 

0 


j 


1,  2 


V. . (t) 
ID 


1 - exp{-hQ [ ( i— 1 ) fQ+t] -H [ ( j-i+1) T -t] } 


(54) 
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(t)  = 1 - exp{ -h  [ (i-1) T +t] -H (kT  -t) -K . . , . } 

13  o o o 3-x^k+l 

k = ] , 2,  (j-i)  (55) 


Kk  = -In  [1  - P(k)] 


(56) 


where  6.  . = 1 if  j-i-2  > 0,  and  6.  . =0  otherwise. 

3-1-2  - j-i-2 

The  probability  of  failure  given  in  Eq.  51  is  for  a single 
aircraft.  Assuming  the  event  of  failure  of  each  aircraft  to  be 
statistically  independent,  the  probability  of  first  failure  in 
a fleet  of  M aircraft  in  [0,  jT  ] interval  is 

PM(j)  = 1 - [1  - P(j)]M  (57) 

2*  Probability  of  Failure  Based  on  Distribution  of  Initial 
Crack  Size  with  No  Cursory  Inspection 

In  this  case,  the  (initial)  crack  of  size  aQ  is  assumed 

to  be  pre-existing  in  the  structure  with  a probability  density 

function  G(aQ).  Although  a number  of  alternative  analytical 

forms  appear  to  be  possible  (Ref.  47),  this  investigation  uses  a 

particular  form  of  G(aQ)  in  the  analysis  which  is  compatible 

with  the  two  parameter  Weibull  density  given  in  Eq . 49. 

The  significance  of  the  compatibility  is  as  follows:  An 

initial  crack  of  size  a increases  its  size  under  loading 

o 

history  and  attains  a specified  size  a (>aQ)  at  time  t . 
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If  the  crack  propagation  law  given  in  Eq.  17  is  used  under 


the  assumption  that  a*  < a^ , then  tQ  and  aQ  are  related  by 


,l-A'/2  l-A'/2 


- a 


t = 
o 


(1  - A' 


/2 ) c ’ (tt/2)  A ,/,2QA 


(53) 


Nov/  interpret  a*  as  the  crack  size  initiated  at  time  t in  the 
crack  initiation  model  considered  in  the  preceding  subsection, 
where  symbol  aQ  is  used  for  the  initiated  crack  size  rather 
than  a*,  however.  Then,  a density  function  G(aQ)  of  initial 
crack  size  aQ  can  be  derived  from  the  density  function  W(tQ) 
of  t with  the  aid  of  the  relationship  given  in  Eq.  58; 


G ( a ) = 
o 


dt 

^W(t  ) 

da  ° 
o 


(59) 


with 


dt 


o 


- a"A'/2 
ao 


da 


o 


* '/2  A ' /b  1-A '/b 
c ' (tt/2)  Q 7 N 

o 


(60) 


If  Eq.  49  is  used  for  W(tQ)  in  Eq.  59,  then  G(aQ)  is  the 
desired  density  function  of  aQ  compatible  with  the  two  para- 
meter Weibull  distribution  of  t . This  compatibility  between 


the  initial  crack  distribution  G(aQ)  and  the  distribution  of 


the  time  to  crack  initiation  W(t  ) is  schematically  shown  in 
Fig. 30  . An  example  of  such  compatible  density  function  is 
plotted  in  Fig.  31  where  the  values  of  Weibull  parameters  as 
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well  as  those  of  the  parameters  involved  in  the  crack  propa- 
gation law  (Eq.  58)  used  for  the  plot  are  indicated.  The 
distribution  function  associated  with  this  density  function  is 
then  plotted  on  the  Frechet  probability  paper  in  Fig. 32  . 

Since  the  Frechet  distribution  is  defined  by  the  density  function 

a*  a -a*-l  -a* 

G(a  ) = (~)  exp  {-(a  / 3* ) ) a 2.  0 (6]) 

° g*  g*  o o 

and  is  not  limited  to  the  right,  the  plot  in  Fig. 32  naturally 

curves  upward  asymptotically  approaching  a vertical  line  at 

aQ  = a*  which  in  this  example  is  equal  to  0.04"  and  is  the  upper 

bound  of  a . 

o 

Under  the  assumptions  that  the  structural  component  has 
a pre-existing  crack  with  a density  function  G(aQ),  that  the 
component  is  subjected  to  a rigorous  inspection  at  the  end  of 
each  period  of  Tq  flight  hours  and  that  the  component  is  replaced 
by  a new  component  if  a crack  is  found  by  the  inspection,  the 
probability  of  failure  of  an  aircraft  within  the  interval 
[0,  jT  ] is  given  by 

oo 

P(j)  = / G(ao)P'  (j)daQ  (62) 


where 
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P'(j)  = 1 - F [a (T) ] exp{ -H (T)  - K^_^} 


j-1  k-1 

- I n F* [a<mT) ]F[a(kT) ]exp{-H(kT)  - K- 
k=2  m=l  3 K 

j-1 

II  F*  [a  (kT)  ] exp{ -H  ( jT)  } (63) 

k=l 

with 

Kj  = - In  [1  - P ( j ) ] (64) 

It  is  noted  that  Eq.  64  involves  P(j)  but  not  P' (j) . This  is 
because  the  distribution  of  initial  crack  size  of  the  replaced 
component  is  assumed  to  be  independent  of  that  of  the  replacing 
component. 

The  comparison  between  Eqs.  62  - 64  and  the  corresponding 
equations  in  the  crack  initiation  model  indicates  that  the 
probability  of  failure  can  be  expressed  with  considerably  less 
analytical  complexities  under  the  assumption  of  initial  crack 
size  distribution,  although  the  integration  in  Eq.  62  must  in 
general  be  carried  out  numerically. 

3.  Probability  of  Failure  Based  on  Distribution  of  Initial 

Crack  Size  with  Cursory  Inspections 

Consider  the  inspection  procedure  where  the  rigorous 
inspection  is  performed  on  the  aircraft  x times  in  its  service 
life  T at  equal  interval  of  T^  with  additional  cursory  inspec- 
tions performed  p times  in  each  rigorous  inspection  interval 
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of  T . Then,  T^  =(p  + 1)Tq,  T = (t  + 1)T1  and  T = (p  + 1)  (x  + 1)TQ 
(see  Fig. 33  ).  Under  these  circumstances,  writing  P(j)  for  the 
probability  of  failure  in  the  interval  [0,  jT-^]  or  [0,  j (p  + 1)TQ], 
one  can  show  that 

oo 

P(j)  = / P'[(p  + l)j]G(a0)dao  (65) 

o 

where 

P'(i)  = 1 " (6  [i-(p+l)n-l]F[a(T0)  ] 

+ [1-6  [ i-  ( p+1)  n-1]  IF  [a(To)]exp[-H(T0)-Ki-]L]  } 
i-1  k-1 

-l  n { 6 [ i — ( p+1) n-m]  F* [a (mTQ) ] 
k=2  m=l 

+ [1—6 [i— ( p+1) n-m] ] F* [a (mTQ) ] } * 

(6[i-(p+l)n-k]F[a(kTo)]  + [1-6 [i-(p+l) n-k] ] 

•F  [a (kT  ) ] }exp{-H(kT  ) - K.  } 

CO  o 1 — K 

i-l 

- n {6 [i- (p+1) n-k] F* [a (kT) ] 
k=l 

+ [1-6  [i-(p  + l)n-k]]Fj[a(kTo)]}exp{-H(iTo)}  (66) 

K.  = -In [ 1 - P(i) ] (67) 

l 
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and 


1 


6 [i  - ( p+1) n - k]  = 


if  there  exists  a positive 
integer  n such  that 
i-k= ( p+1) n (68) 

0 otherwise 


4.  Probability  of  Failure  of  Proof  Tested  Aircraft 

To  include  the  effect  of  the  proof  load  test,  define  a 
probability  density  function  G' (aQ)  for  the  initial  crack 
distribution  in  the  following  fashion  (Refs.  48  and  49) 


(G(ao)/  /o°P G(ao)dac 


G' (ao}  ~ 


' 0 


0 < a < a 

o — Op 


(69) 


a < a 
op  o 


where  a is  the  truncation  level  on  G(a  ) as  a result  of  the 
op  o' 

proof  test.  The  probability  of  failure  can  then  be  determined 
by  the  same  procedures  as  described  in  8.2  or  8.3  depending  on 
whether  or  not  the  cursory  inspections  are  performed.  The 
necessary  modifications  consist  of  replacing  G(aQ)  by  G'(aQ)  and 
the  infinite  upper  bound  by  a in  the  intervals  in  Eqs.  62  and  65. 

It  appears  difficult  at  this  time  to  consider  the  effect  of 
the  proof  load  test  within  the  framework  of  the  crack  initiation 
model  described  in  8.1.  Furthermore,  it  is  pointed  out  that  Eq . 
57  also  indicates  the  probability  of  first  failure  in  a fleet 
of  M aircraft  in  the  interval  [0,  j T ] for  the  cases  of  8.2  and 
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in  the  interval  [0,  jT^]  in  the  case  of  8.3.  In  the  above, 
8.1,  8.2  etc.  indicate  respectively  subsections  1,  2,  etc,  of 
this  section  (Section  VIII) . The  same  notation  is  used  for 
other  subsections  throughout. 
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IX.  Sensitivity  Analysis 
1.  Model  Sensitivity 

In  the  preceding  sections,  a number  of  analytical  models 
are  introduced  with  respect  to  random  stress  processes,  crack 
propagation  laws,  rise  and  fall  statistics  of  the  stress  process, 
residual  strengths,  expected  failure  rate  and  inspection  pro- 
cedures for  the  ultimate  purpose  of  evaluating  the  probability 
of  aircraft  failure.  Obviously,  the  resulting  probability  of 
failure  depends  on  these  analytical  models.  Since,  however,  they 
are  contructed  on  the  basis  of  the  current  (and  therefore  not 
necessarily  perfect)  state  of  the  engineering  knowledge  on  the 
subjects  involved  with  the  due  consideration  for  analytical  and 
numerical  tractability , it  is  highly  desirable  to  investigate 
the  sensitivity  of  the  probability  of  failure  to  these  models. 

Such  sensitivity  studies  could  be  performed  in  terms  of  numerical 
examples  for  those  items  for  which  competing  or  alternative 
models  are  proposed.  For  example,  significant  difference  in  the 
probability  of  failure  may  result  depending  on  whether  one  choses 

(a)  a single  Gaussian  or  a composite  Gaussian  for  a com- 
ponent process  (such  as  stress  process  due  to  clear  air 
turbulence) , 

(b)  a Gaussian  process  or  a non-Gaussian  process  for  maneuver, 

(c)  counting  method  I or  method  II  in  conjunction  with  rise 
and  fall  statistics  (see  Fig.  17), 
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(d)  slow  crack  growth  design  or  fail-safe  design, 

(e)  clumpsize  equal  to  one  or  any  other  value, 

(f)  crack-initiation  model  or  pre-existing  crack  model, 

(g)  controlled  maneuver  or  uncontrolled  maneuver, 

(h)  perform  inspection  or  do  not  perform  inspection 

(i)  perform  proof  load  test  or  do  not  perform  proof  load 
test . 

Although  the  computer  programs  are  written  so  that  all  these 
items  except  for  item  (e)  can  be  dealt  with,  the  sensitivity 
studies  are  actually  performed  only  on  items  (d) , (f) , (g) , (h) 

and  (i)  and  the  results  are  summarized  in  terms  of  graphical 
comparisons  (see  Sec.  XI).  Such  results  are  highly  useful  in 
numerically  assessing  the  significance  of  various  models  des- 
cribed above. 

2.  Parameter  Sensitivity 

After  the  decision  is  made  to  choose  a particular  model 
among  the  alternatives  proposed  for  each  of  the  items  indicated 
above,  it  is  not  unusual  to  face  the  uncertainties  of  statistical 
and  other  origins  under  which  the  values  are  to  be  assigned  to 
the  parameters  of  the  model.  Then,  the  sensitivity  studies  are 
again  desirable  with  respect  to  these  parameters  for  the  purpose 
of  identifying  more  important  parameters  in  terms  of  their 
contributions  to  the  probability  of  failure. 

Writing  these  parameters  as  Xw  X0,  X , the  probability 

\ i n 
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of  failure  P can  formally  be  written  as  a function  of  X^ , X^, 


, X^  once  the  analytical  models  to  be  used  are  decided  upon 


P = P(X.  , X , X ) 

12  n 


(70) 


In  reality,  the  dependence  of  P on  the  parameters  is  obviously 
very  complex  as  seen  from  the  analysis  developed  in  the  pre- 
ceding sections. 

Treating  the  uncertainties  associated  with  these  parameters 
as  if  they  were  all  of  statistical  origin  (or  treating  X^,  X^,  . 
X^  as  random  variables) , the  sensitivity  of  the  probability  of 
failure  to  each  of  these  parameters  is  evaluated  with  the  first 
order  and  second  moment  approach.  Within  the  framework  of  this 
approach,  the  most  reliable  value  of  a parameter  X^  is  treated 
as  its  expected  value  and  written  as  X^ , whereas  its  uncertainty 
is  expressed  in  terms  of  coefficient  of  variation  Vx  = oYi/  x\ 


'X, 


with  cv  being  standard  deviation  of  X..  The  expected  value 

1 

P and  coefficient  of  variation  Vp  of  P are  then  obtained  from 


P = P(xn  , X . X ) 

12  n 


(71) 


V2  = T a2V2 
P i=l  i Xi 


(72) 


where  is  referred  to  as  sensitivity  index  and  is  given  by 
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a . 
1 


(73) 


1 


where  partial  derivatives  are  evaluated  at  X,  = X^,  (i  = 1,  2,..., 
n) . If  P and  Vp  are  interpreted  as  representing  the  most  reliable 
value  of  P and  its  uncertainty  respectively,  the  sensitivity 
index  a_^  indicates  the  contribution  of  the  uncertainty  associ- 
ated with  parameter  to  the  uncertainty  of  P through  Eq.  72. 

Dealing  with  a transport  aircraft,  twenty  five  (25)  para- 
meters are  investigated  for  the  sensitivity  analysis  and  are 
listed  in  Table  - 5;  Table  3 also  shows  their  most  reliable 
values  as  expected  values  and  corresponding  sensitivity  indices 
squared  (a?)  for  crack  initiation  model  under  slow  crack  growth 
design.  Table  4 for  crack  initiation  model  under  fail  safe  design 
and  Table  5 for  pre-existing  crack  model  under  fail  safe  design. 
These  results  are  all  obtained  under  the  assumption  of  five 
equally  spaced  rigorous  inspections  performed  during  the  service 
life  of  15,000  hours.  As  mentioned  earlier,  the  dependence  of 
P on  Xp  X , X^  is  highly  complex  and  analytical  evalua- 

tions of  partial  derivatives  in  Eq.  73  are  in  general  not 
possible.  Therefore,  these  derivatives  are  numerically  evaluated 
in  principle  as 


={p(x1 


, X . + Ax  • . . . ,X  )-P(X.  ,xo 
i i n 12 

(i  = 1,  2, 


,X. , . . . ,X  ) }/AX. 
x n l 

n)  (74) 
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The  procedure  involved  in  such  numerical  evaluation  is  illust- 
rated in  Figure  34  dealing  with  aQ,  the  first  parameter  in 
Table  2 as  an  example. 
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X. 


Significance  of  Full-Scale  Test 


Fatigue  testing  under  loading  sequences  representative  of 
expected  operational  conditions  obviously  provides  realistic 
life  estimates.  These  estimates  are,  however,  usually  obtained 
from  full-scale  tests  involving  a single  or  at  most  two  proto- 
types, since  the  cost  is  prohibitive  even  for  a small  number 
of  replications.  This  subsection  therefore  attempts  to  inves- 
tigate the  significance  of  such  full-scale  testing  from  the 
view  point  of  aircraft  reliability  estimation.  In  this  respect, 
the  analysis  recently  developed  in  Ref.  50  is  closely  followed 
where  the  time  to  first  failure  t-^ , is  related  to  a reliability 
figure  on  the  basis  of  the  result  of  full-scale  testing  without, 
however,  considering  the  effect  of  inspections. 

Although  the  full-scale  testing  is  performed  usually  only 

on  a single  or  at  most  two  prototypes  as  just  mentioned, 

consider  fo  r the  purpose  of  theoretical  development  that  n 

such  prototypes  are  tested  with  observed  lives  (times  to  failure) 

t , , t t . These  observations  constitute  a sample  of 

ol  o2  on  ^ 

size  n taken  from  the  density  function  W(tQ)  of  the  time 
to  failure  t of  individual  aircraft  (or  structural  components) . 
Since,  as  will  be  seen  below,  this  sample  is  needed  to  estimate 
a measure  of  location  of  the  distribution  of  t , any  density 
function  that  provides  a reasonable  fit  to  this  sample  at 

a central  range  of  distribution  may  be  assumed  for  W(tQ) . For 
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convenience,  however,  the  two  parameter  Weibull  distribution 
with  the  density  function  given  in  Eq.  49  is  assumed  for  tQ. 
The  time  to  first  failure  t^  in  a fleet  of  size  m is  then 
interpreted  as  the  smallest  among  the  sample  of  size  m taken 
from  the  Weibull  distribution  with  the  density  W ( tQ)  given  in 
Eq.  49.  Hence,  the  distribution  function  of  t^  is  in  general 
one  of  the  extreme  value  distributions  and  in  this  case  it  is 
again  a Weibull  by  virtue  of  stability  of  the  extreme  value 
distribution.  In  fact,  it  can  be  shown  to  be 


W 


- ( ) exp{- ( t /B , ) } 

&1  3i 


(75) 


where,  however,  the  scale  parameter  3^  is  given  by  3^  = $/m 
while  the  shape  parameter  remains  to  be  a (a  and  3 are  shape  and 
scale  parameter  of  the  Weibull  density  for  t ) . 

Define  now  the  scatter  factor  S as 


S = 


B/t1 


(76) 


where  3 is  the  maximum  likelihood  point  estimator  of  the  scale 
parameter  3 from  a sample  of  size  n (tQ^,  t - , ....,  t ). 

The  point  estimator  3 for  n observations  t is  given  by 


6 


1 

[— 

n 


n 


I t“  ]1/a 

i=l  oi 


(77) 
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provided  a is  known.  The  distribution  function  f^(S)  can  then 
be  shown  to  be 


f x ( $)  d$  = 


n 

n a 


f(n)  3 


rB',  an-1  a.  ~ 

• (p  exp[~n(p)  ]d3 


(78) 


Using  Eqs.  75  and  78,  one  obtains  the  density  function  of  the 

✓\ 

quotient  S = $/t  as 


oo 

f2(S)  = /0fl(Stl} ‘WV^l 


(79) 


which  after  some  manipulation  becomes 


n+lnan-l 
amn  S 


f2(S)  = 


(m  + nSa)n+1 


(80) 


By  integration,  the  distribution  function  of  S is  obtained  as 

a 

® n 

F_(S)  = [ ~]  (81a) 

* (m/n)  + Sa 

and  this  probability  is  identical  to  the  reliability  R that  the 
time  to  first  failure  t,  will  be  at  least  equal  to  the  maximum 
likelihood  estimate  defined  by  Eq.  77  divided  by  a scatter 
factor  S; 

Sa 

R = [ cr]n  (81b) 

(m/n)  + S 
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The  values  of  a representative  for  different  structural 
materials  and  for  aircraft  of  fighter  and  transport  type  are 
summarized  in  Ref.  50  on  the  basis  of  existing  data  and  listed 
here  in  Table  6.  The  values  of  the  scatter  factor  associated 
with  reliability  levels  R = 0.5,  0.9  and  0.99  are  evaluated 
from  Eq . 81b  for  more  conservative  values  of  a = 2,  3,  4 and  5. 
The  results  are  shown  in  Tables  7-9  for  m = 3,  25,  100,  250 
and  1,000  and  m = 1 and  3.  It  is  important  to  observe  from 
Tables  7-9  that  practically  no  change  materializes  by 
increasing  n from  1 to  3 . This  fact  is  highly  significant,  as 
emphasized  in  Ref.  50,  since  it  means  that  the  full-scale 
testing  performed  on  two  or  three  prototypes  offers  practically 
no  advantage  over  the  testing  of  a single  prototype  as  far  as 
the  reliability  estimation  with  respect  to  the  time  to  first 
failure  is  concerned.  It  is  also  important  to  realize  that  the 
values  of  the  scatter  factor  are  within  a practical  range  even 
at  a reliability  level  of  .99  particularly  when  the  fleet  size 
is  not  too  large.  If  the  minimum  life  T*  is  introduced  in  the 
Weibull  distribution  as  indicated  in  Eq.  50,  it  can  be  shown 
(Ref.  50)  that  the  values  of  the  scatter  factor  will  be  reduced 
significantly  even  for  T*  as  small  as  0.058. 

The  results  described  above  provide  the  theoretical  bases 
on  which  a realistic  method  of  reliability  demonstration  can  be 
established  as  an  enforcable  part  of  the  certification  and 
procurement  procedure.  Even  more  significant  is  the  potential 
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of  this  approach  in  evaluating  the  possible  amount  of  reduc- 
tion in  the  scatter  factor  when  the  aircraft  are  subjected  to 


the  periodic  inspection  as  described  in  Sec.  VII, 
approach  must  in  essence  use  the  crack  initiation 
of  its  mathematical  development  using  the  Weibull 


although  the 
model  because 
distribution. 
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XI.  Numerical  Examples  and  Discussion 


Numerical  examples  are  carried  out  for  long-life  aircraft 
such  as  transports  and  bombers  using  the  crack  initiation  model 
under  slow  crack  growth  and  fail  safe  design  and  also  using  the 
pre-existing  crack  model  under  fail  safe  design.  With  the 
crack  initiation  model,  a two-parameter  Weibull  density  is  used 
for  the  time  to  crack  initiation  while  the  corresponding  (com- 
patible) distribution  is  used  with  the  pre-existing  crack  model. 
Tables  3-5  show  the  parameter  values  used  for  the  numerical 
computation  in  the  column  under  "Mean  Value."  Some  of  these 
values  are  assigned  on  the  basis  of  experimental  evidence, 
others  are  chosen  to  be  consistent  with  conventional  static 
design  practice  and  still  others  are  on  the  basis  of  engineering 
judgement.  For  example,  the  values  of  material  constants  are 
consistent  with  those  of  structural  aluminum,  the  parameters 
associated  with  the  detectable  crack  size  are  within  the  range 
found  in  available  experiments,  and  the  power  b in  the  crack 
propagation  law  is  assumed  to  be  4.0  because  of  the  availability 
of  the  experimental  data  under  random  loading  only  for  b = 4. 

The  parameters  involved  in  the  formulation  of  stress 
processes  are  more  difficult  to  estimate  since  not  much  data 
are  available  in  terms  of  stress  history.  The  usual  load  data 
in  terms  of  g cannot  directly  be  used  without  specific  knowledge  Q.f 
the  structural  system  under  consideration.  The  coefficient  A 
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1 


appearing  in  the  rise  and  fall  statistics  is  obtained  from 
the  result  of  simulation  with  b = 4 and  counting  method  II 
(Fig.  17);  A = A14  = A^  = 115  for  S c(t)  (CAT)  and  SgT(t) 
(thunderstorm)  assuming  that  they  have  narrow-band  characteris- 
tics equivalent  to  those  simulated  processes  with  Bc  = 0.6  - 0.8 
(see  Fig.  14  and  19),  while  A = A^  = 80  for  Sm(t)  (maneuver) 
from  the  solid  curve  in  Fig.  25  assuming  that  symmetric  excee- 
dance curves  apply  to  transports  (Fig.  2).  The  number  of  stress 
cycles  per  hour  such  as  N and  represents  the  expected  rate 
of  zero  upcrossing  since  the  second  counting  method  is  used  in 
this  study.  Therefore,  the  evaluation  of  these  values  is  obviously 
possible  when  the  spectral  density  functions  of  the  underlying 
stress  processes  are  specified.  Because  of  the  fact  that  the 
spectral  densities  of  "stress"  processes  are  usually  not  well 
specified,  Ng  and  are  treated  as  independent  parameters.  It 
is  emphasized  that  this  does  not  invalidate  either  the  inter- 


pretation or  the  result  of  the  simulation  to  evaluate  the 
coefficient  A since  one  can  always  adjust  the  time  scale  so 
that  a desired  value  of  Ng  or  Nm  is  obtained  without  altering 
the  value  of  A. 

The  choice  of  p^  = 0.495,  p2  = 0.005  and  p^  = 0.5  is  made 
because  such  a choice  appears  to  represent  one  of  realistic 


loading  conditions.  Also,  the  values  of  a^,  Op2 


and  cr  _ as 


indicated  in  the  tables  are  judged  to  be  reasonable. 


Under  the  loading  and  other  conditions  corresponding  to 
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these  parameter  values  specified  in  Table  4,  a crack  of  size 
0.04",  once  initiated,  grows  under  fail  safe  design  as  shown 
in  Fig.  35  where  the  corresponding  residual  strength  R(tn)  and 
failure  rate  h(t^)  are  also  plotted.  Note  that  the  crack  be- 
comes unstable  after  4,000  - 5,000  flight  hours. 

Figs.  36  - 39  illustrate  the  effect  of  some  of  the  model 
sensitivity  mentioned  in  Sec.  IX.  The  probabilities  of  failure 
are  all  in  terms  of  the  probability  of  first  failure  in  a fleet 
of  size  50  for  service  life  of  15,000  flight  hours.  Therefore, 
the  probability  of  failure  of  any  aircraft  is  in  approximation 
one  fiftieth  of  the  value  indicated. 

Figs. 36  (a)  and  (c)  show  the  difference  between  the  pro- 

bability of  first  failure  resulting  from  the  pre-existing  model 
and  that  resulting  from  the  crack  initiation  model  both  under 
fail  safe  design.  The  pre-existing  model  and  the  crack  initia- 
tion model  are  labeled  as  G(aQ)  method  and  W(tG)  method  for 
simplicity.  Figs.  36  (b)  and  (d)  show  the  same  difference  under 
the  slow  crack  growth  model.  When  no  inspection  is  performed, 
the  pre-existing  crack  model  and  the  crack  initiation  model 
produce  practically  identical  probability  of  first  failure  under 
both  fail  safe  and  slow  crack  growth  designs.  The  reason  is  as 
follows;  when  G(a0)  and  W(tQ)  are  compatible,  the  reduction  in 
the  ultimate  strength  during  the  time  period  t in  the  pre- 
existing model  is  negligibly  small,  at  least  for  the  numerical 
examples  considered,  compared  with  the  ultimate  strength  RQ 
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which  remains  constant  until  the  end  of  the  same  time  period 
under  the  crack  initiation  model.  Fig.  36  also  illustrates 
the  effect  of  rigorous  inspection:  In  all  cases,  the  proba- 

bility of  first  failure  decreases  as  the  number  N of  rigorous 
inspections  increases.  However,  no  further  improvement  is 
expected  for  N larger  than  14.  So  long  as  N < 14,  the  improve- 
ments achieved  by  the  same  number  of  rigorous  inspections  are 
more  significant  for  the  pre-existing  crack  model  than  for  the 
crack  initiation  model  under  both  fail  safe  and  slow  crack 
growth  design.  This  is  due  to  the  fact  that,  as  the  pre-existing 
crack  aQ  grows,  it  may  be  detected  possibly  at  the  early  stage 
of  growth  with  a size  much  smaller  than  a*  (=  0.04"  in  this 
example)  and  the  cracked  structural  component  replaced  thus  pro- 
viding better  chances  for  containing  the  crack  size  less  than 
a*  throughout  the  life,  while  in  the  crack  initiation  model  the 
cracks  are  found  only  after  they  are  initiated  with  the  size  of 
a*.  Comparison  between  Figs.  36  (a)  and  (b)  and  between  (c) 
and  (d)  show  that  the  fail  safe  design  produces  safer  aircraft 
with  both  pre-existing  crack  model  and  crack  initiation  model 
at  least  for  this  numerical  example.  It  is  important  to  note 
that  all  the  curves  in  Fig.  36,  in  particular,  the  curve  with 
no  inspection,  show  faster  rates  of  increases  in  probability 
of  failure  after  about  4,000  - 5,000  flight  hours  and  that  this 
is  the  range  at  which  a(tn)  start  to  increase  significantly 
as  mentioned  previously  with  respect  to  Fig.  35.  The  same 
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trend  is  observed  in  Figs.  37  - 41  that  will  follow.  The 
exceptions  are  those  under  the  favorable  effect  of  proof  load 
test  in  Fig.  39. 

Figs.  37  and  38  indicate  the  effect  of  cursory  inspections 
under  the  various  assumptions  of  the  crack  size  a^ 1 that  can 
be  detected  by  them.  These  results  are  for  the  pre-existing 
crack  model  under  fail  safe  design.  Fig.  37  shows  that  even 
when  a^ ' is  as  small  as  0.1",  the  improvements  on  the  reliabi- 
lity achieved  by  the  cursory  inspections  are  insignificant  as 
long  as  the  number  of  cursory  inspections  is  small.  Each  dia- 
gram represents  a combination  of  different  numbers  of  rigorous 
and  cursory  inspections.  As  shown  in  Fig.  38,  however,  signifi- 
cant improvements  can  be  achieved  by  increasing  the  number  of 
cursory  inspections  considerably.  Even  then,  however,  the 
detectable  size  must  be  as  small  as  0.1". 

The  results  given  in  Fig.  39  is  significant  since  they 
indicate  generally  favorable  effect  of  proof  load  tests  (under 
both  fail  safe  and  slow  crack  growth  designs  with  the  pre-exist- 
ing model)  particularly  when  a reasonable  number  of  rigorous 
inspections  are  performed;  for  example,  5 inspections  in  Figures 
39  (b)  and  (d) . It  should  be  noted  that  there  apprears  to  be 

a critical  truncation  value  aQp  such  that  proof  load  tests 
producing  a truncation  less  than  such  a critical  value  would 
significantly  improve  the  reliability.  Campare , for  example, 
the  probability  of  first  failure  for  15,000  flight  hours  for 
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a = 0.030  with  that  for  a = 0.025  in  both  Figs.  30  (a) 
op  op 

and  (c)  . 

The  parametric  sensitivity  studies  are  performed  on  all 
twenty- five  parameters  listed  in  Tables  3-5  using  the  method 
described  in  subsection  9.2  and  following  the  procedures  illus- 
trated in  Fig.  34.  In  fact.  Fig.  40  (d)  shows  how  the  basic 
probability  values  made  use  of  in  Fig.  34  for  the  computation 

of  sensitivity  index  a are  obtained.  Other  diagrams  in  Fig. 

i 

40  show  corresponding  probability  values  for  the  parameters 

(AK^) and  g (stress  equivalent  of  one  g load)  both  for  Table 
TH 

5 and  a for  Table  4.  The  resulting  values  of  a2  are  listed 
o i 

under  the  column  a2  in  Tables  3-5  (all  for  the  case  of  five 

i 

rigorous  inspections).  In  all  cases,  the  parameter  g is  most 

b 1/b 

important,  followed  by  (AK  ) . Other  more  significant 

TH  — 

4 

parameters  are  found  to  be  X ' , a_  , a or  a*,  Z , etc.  For 

Ci  o 

fail  safe  design,  y and  £ are  also  important. 

o 

Fig.  41  indicates  the  probability  of  first  failure  when  no 

~5  1/b 

threshold  is  considered  for  (AK  ) with  log(da/dn)  being 

b 1/b 

related  to  log(AK  ) through  a straight  line  (CD  in  Fig.  8). 

Figs.  42  and  43  give  the  overall  effect  of  the  number 
of  inspections  for  the  pre-existing  crack  model  and  the  crack 
initiation  model  respectively. 

For  fighters,  numerical  examples  are  given  with  the 
crack  initiation  model.  In  one  case,  structural  details  made 
of  steel  are  considered  under  slow  crack  growth  design. 
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The  parameters  listed  in  Table  3 are  used  for  numerical  computa- 
tion except  for  g = 15.5  ksi,  P = 1.0,  A - 12,  y = 210  ksi, 
— b 1/b  _ 3 _9  34  o 

(AK  ) =7.0  ksi/in  and  c = 3 x 10  . Also,  in  this  case 

TH 

asymmetric  processes  constructed  from  single  Gaussian  processes 

as  introduced  in  Eq.  30  are  used  for  the  stress  processes  Sm(t) 

for  maneuver  with  r = 0.3  (this  is  why  A = 12,  see  Fig.  25), 

34 

o=o=  between  1.55  g and  0.5  g.  The  value  a = 1.55  g is 
S p P 

consistent  with  the  exceedance  curves  given  in  Fig.  3 and  appears 

to  represent  the  upper  bound  for  a , while  the  value  o = 0.5  g 

P P 

is  judged  to  represent  less  strenuous  maneuver  conditions. 

The  results  are  shown  in  Figs.  44  and  45.  Fig.  44  shows  the 

effect  of  the  number  of  rigorous  inspections  while  Fig.  45 

the  effect  of  a . Figs.  46  - 48  indicate  similar  results 
S 

for  structural  components  made  of  aluminum.  The  parameters 
are  identical  to  those  listed  in  Table  3 and  4,  depending  on 
whether  slow  crack  growth  design  or  fail  safe  design  is  con- 
sidered. The  exceptions  are  g = 4.2  ksi,  P =1.0  and  A = 

3 34 

12.  The  same  asymmetric  stress  processes  are  used;  r = 0.3 

(and  therefore  A = 12)  and  o = a = 0.5  g - 1.55  g.  Fig. 

34  S p 

46  shows  the  effect  of  the  number  of  rigorous  inspections  under 
fail  safe  design  while  Figs.  47  and  48  indicate  the  effect 


of  a under  fail  safe  design  and  slow  crack  growth  design 
S 

respectively.  It  is  pointed  out,  however,  that  the  results 


indicated  in  Figs.  45,  47  and  48  do  not  really  represent  the 


effect  of  a since  the  identical  Weibull  density  has  been  used 
S 

for  the  time  t to  crack  initiation  irrespective  of  a values 
o S 

while  in  reality  the  time  t obviously  depends  on  o . 

o S 
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XII.  Conclusion  and  Recommendation  for  Further  Study 

The  probability  of  aircraft  failure  P^  may  be  considered 
as  function  of  a large  number  of  variables.  For  example, 

Pf  = Pf (M , F,  D,  G,  S,  T,  I,  P,  A) 

where 

M = Material  selection 

F = Fabrication 

D = Design  practice 

G = Geometric  configuration 

S = Mission  spectra 

T = Full  scale  and  other  testing 

I = Inspection  procedures 

P = Proof  test 

A = Analysis  methods 

In  the  present  study,  a major  effort  has  been  made  to 
incorporate  into  the  reliability  evaluation  scheme,  the  effect 
of  material  selection,  geometrical  configuration  (fail  safe  or 
slow  crack),  mission  spectra,  full  scale  testing,  inspection  pro- 
cedures, and  proof  load  test.  The  effect  of  design  practice 
implicitly  appears,  for  example,  through  the  stress  equivalent 
of  one  g load  as  the  design  usually  equates  the  ultimate  strength 
to  the  limit  load  in  terms  of  g multiplied  by  a safety  factor. 
Also,  the  analysis  method,  which  determines  the  accuracy  of  the 
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(calculated)  stress  equivalent  of  one  g load  as  well  as  of 
(calclulated)  strsss  processes,  influences  the  reliability. 

Although  their  roles  are  implicit  in  the  analysis  presented  here, 
these  are  important  items  from  the  view  point  of  reliability 
evaluation  since  the  probability  of  failure  has  been  found 
to  be  extremely  sensitive  to  the  parameter  g . 

Although  this  study  has  provided,  it  is  hoped,  an  adequate 
analytical  framework  for  the  current  effort  toward  implementation 
of  the  reliability-based  design  criteria  for  USAF  aircraft, 
there  are  a number  of  other  items  on  which  further  studies 
are  recommended.  They  include  the  establishment  of 

(1)  effect  of  inspection  procedures  on  the  scatter  factor 

(2)  more  reliable  crack  propagation  law,  particularly  values  of 
(AKd)th  and  X',  under  random  loading 

(3)  more  reliable  probability  values  for  detecting  cracks 
of  various  size  by  NDI  techniques 

(4)  more  realistic  model  for  residual  strength  for  both 
slow  crack  growth  and  fail  safe  design 

(5)  range  of  the  parameter  values  associated  with  stress 
processes  from  existing  data  depending  on  loading  conditions 
and  structural  details,  particularly  for  fighters  under 
maneuver  loading 

(6)  more  explicit  dependence  of  the  distribution  of  time 
to  failure  on  the  mission  spectra  and  the  crack  size  at 
the  time  of  crack  initiation  with  respect  to  the  crack 


66 


initiation  model 


(7)  more  reliable  distribution  functions  for  initial 
crack  size  with  respect  to  the  pre-existing  crack  model 

(8)  effect  of  having  a large  number  of  critical  locations 
on  the  reliability  value  of  an  aircraft 

(9)  optimum  number  of  periodic  inspections  on  the  basis 
of  cost-effectiveness  consideration. 

Although  all  these  items  are  important,  the  items  (1)  - (5),  when 
accomplished,  will  have  immediate  impact  on  the  implementation 
of  the  reliability-based  design  and  procurement  procedures. 
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Figure  5.  Fatigue  Crack  Growth  Rates  at  Low  Stress  Intensities 
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Figure  6.  Crack  Growth  Rates  under  Random  Loading  (Ref.  36) 
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Crack  Growth  Rates  under  Sinusoidal  and  Random  Loading 
(Carbon  Steel)  (Ref.  38) 


Figure  8.  Crack  Growth  Rate  Model  Used  in  Analysis 
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Figure  16.  Sample  History  ($c“  0.9) 
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Figure  17.  Counting  Methods 
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Figure  23.  Sample  History  (Maneuver;  Asymmetric) 


Figure  25.  Coefficient  A as  a Function  of  Power  b (Maneuver) 
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Figure  30.  Compatibility 
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Figure  35.  Residual  Strength,  Crack  Size  and  Failure  Rate  as  Functions  of  tn 
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Figure  36.  Effect  of  Rigorous  Inspections 
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Figure  39.  Effect  of  Proof  Load  Test 
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Figure  40.  Sensitivity  Study  (Examples) 
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Figure  42.  Effect  of  Number  of  Inspections  (W(tQ)  Method) 
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PROBABILITY  OF  FAILURE  AT  15000  HOURS 


Figure  43.  Effect  of  Number  of  Inspections  (G(aQ)  Method) 
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Figure  44.  Probability  of  First  Failure  of 

a Fleet  of  50  Fighters;  Effect  of 
Number  of  Inspections  (Steel) 
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Figure  45,  Probability  of  First  Failure  of  a Fleet 
of  50  Fighters;  Effect  of  Go  (Steel) 
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TABLE  1.  IRREGULARITY  RATIO  R 

r = — j (l-e3)/  J ( i-6)  (i-e5) 


B 

Theory 

Simulation 

Error (%) 

0 

.7454 

.7503 

.65 

25 

.8476 

. 8485 

.10 

50 

. 9371 

. 9364 

. 07 

75 

. 9868 

. 9884 

.16 

80 

.9919 

.9898 

.21 

90 

. 9982 

. 9984 

.02 
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TABLE  2 


EXPECTED  VALUE  OF  b-TH  POWER  OF  RISE  AND  FALL 


(a)  b = 1 


3c 

Simulation 
Ref.  29 (Exact) 
Ref.  30 


0 0.25 
1.87  2.12 
1.87  2.11 


0.5  0.75 
2.35  2.46 
2.35  2.48 


0.8  0. 

2.47  2. 


(b)  b = 4 


3c 

0 

0.25 

0.5 

0.75 

0.8 

0.9 

Simulation 

48 . 2 

77 . 0 

108 

121 

125 

126 

Ref.  29 

43.6 

63.1 

81.3 

133 

- 

- 

Ref.  30 

- 

.. 

... 

- 

127 

128 
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L n to 


TABLE  3.  SENSITIVITY  ANALYSIS  (SLOW  CRACK  GROWTH) 


Parameter 

Mean  Value 

2 , r 
a.  (N=5 

X 

1. 

a : Crack  size  initiated  at  t 
o o 

0.04  in 

50.21 

2. 

c:  Material  constant 

3.0  x 10-7 

22.96 

3. 

Ng : Number  of  gust  load  cycles  per  hour 

600 

0.069 

4. 

N : Number  of  maneuver  load  cycles 
m per  hour 

60 

0.156 

5. 

N : Number  of  G-A-G  load  cycles  per  hour 
z 

0.5 

4.866 

6. 

A = A 

14  24 

115 

8.216 

7. 

A34 

80 

1.602 

8. 

P^:  Fraction  of  flight  hours  in  maneuver 

0.5 

0.022 

9. 

0 : Intensity  of  CAT 

0.07g  ksi 

111.89 

10. 

a^:  Intensity  of  thunderstorm 

0.18g  ksi 

3.887 

11. 

aC3:  Intens^-t^  of  maneuver 

0.10g  ksi 

12.75 

12. 

? : Average  of  fourth  power  of  G-A-G 
cycles  Z 

(1.5  g)4 

102.35 

13. 

g:  Gravitational  acceleration 

10  ksi 

509.98 

• 

i — 1 

K : Critical  value  of  stress  intensity 
c factor 

60  ksi/in 

0.0750 

15. 

V : COV  of  ultimate  strength 
o 

0.056 

0.00034 

16. 

y : Ultimate  strength 
o 

57  ksi 

0.00112 

17. 

a:  Shape  parameter 

4.0 

18.84 

18. 

3:  Scale  parameter 

30,000  hr  s 

15.89 

19. 

U-^:  Probability  of  inspecting  cracked 
details 

1.0 

60.99 

20. 

a^:  Maximum  undetectable  crack  size 

0.02  in 

0.251 

21. 

a^:  Minimum  crack  definitely  detectable 

0.3  in 

0.833 

22. 

m:  Parameter  appearing  in  expression 
for  U2  (a) 

0.125 

2.973 

23. 

(AK^)i^:  Threshold  value 

in 

1.5  ksi/in 

240.79 

24. 

X*:  Slope  of  crack  propagation  law 

10.0 

127.88 

25. 

X:  Slope  of  crack  propagation  law 

2.5 

41.80 
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TABLE  4.  SENSITIVITY  ANALYSIS  (FAIL  SAFE) 


Parameter 

Mean  Value 

a?  (N= 

1. 

a : Crack  size  initiated  at  t 
o o 

0.04  in 

76.69 

2. 

c:  Material  constant 

3.0  x 10_7 

11.38 

3. 

N : Number  of  gust  load  cycles  per  hour 
g 

600 

1.291 

4. 

N : Number  of  maneuver  load  cycles 
m per  hour 

60 

0.252 

5. 

N : Number  of  G-A-G  load  cycles  per 
z hour 

0.5 

5.257 

6. 

A = A 

14  24 

115 

7.409 

7. 

A34 

80 

0.448 

8. 

P : Fraction  of  flight  hours  in  maneuver 
3 

0.5 

0.443 

9. 

as  Intensity  of  CAT 

0 . 07g  ksi 

151.84 

10. 

a^2:  Intensity  of  thunderstorm 

0 . 18g  ksi 

49.42 

11. 

a : Intensity  of  maneuver 
C3 

0.10g  ksi 

A 

10.21 

12. 

Z^ : Average  of  fourth  power  of  G-A-G 
cycles  Z 

(1.5  g) 

86.85 

13. 

g:  Gravitational  acceleration 

10  ksi 

505.58 

14. 

Residual  strength  ratio 

0.43 

95.96 

15. 

V : COV  of  ultimate  strength 
o 

0.056 

0.220 

16. 

y : Ultimate  strength 
o 

57  ksi 

98.97 

17. 

a:  Shape  parameter 

4.0 

19.17 

18. 

g:  Scale  parameter 

30,000  hrs 

15.95 

19. 

U : Probability  of  inspecting  cracked 
details 

1.0 

101.18 

20. 

a^:  Maximum  undetectable  crack  size 

0.02  in 

0.224 

21. 

a : Mimimum  crack  size  definitely 
2 detectable 

0.3  in 

1.060 

22. 

m:  Parameter  appearing  in  expression 
for  (a) 

0.125 

3.151 

23. 

(Ak'3)^*3:  Threshold  value 
1 H 

1.5  ksi/in 

227.70 

24. 

X':  Slope  of  crack  propagation  law 

10.0 

52.54 

25. 

X:  Slope  of  crack  propagation  law 

2.5 

35.01 
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TABLE  5.  SENSITIVITY  ANALYSIS  (FAIL  SAFE) 


Parameter 

Mean  Value 

a?  (N=5) 

1. 

a*:  Upper  bound  of  G(aQ) 

0.04  in 

129.13 

2. 

c:  Material  constant 

3.0  x 10~7 

20.40 

3. 

N^:  Number  of  gust  load  cycles  per  hour 

600 

1.208 

4. 

N : Number  of  maneuver  load  cycles  per  hour 
m 

60 

0.420 

5. 

N : Number  of  G-A-G  load  cycles  per  hour 
z 

0.5 

15.26 

6. 

A14  = A24 

115 

15.89 

7. 

A3  4 

80 

0.911 

8. 

P : Fraction  of  flight  hours  in  maneuver 
3 

0.5 

0.098 

9. 

ac^:  Intensity  of  CAT 

0.07g  ksi 

222.45 

10. 

a : Intensity  of  thunderstorm 

0.18g  ksi 

81.00 

11. 

: Intensity  of  maneuver 

0.10g  ksi 

21.18 

12. 

Z^:  Average  of  fourth  power  of  G-A-G 
cycles  Z 

(1.5g) 4 

180.22 

13. 

g:  Gravitational  acceleration 

10  ksi 

1322.31 

14  . 

£ : Residual  strength  ratio 

0.43 

93.13 

15. 

Vq:  COV  of  ultimate  strength 

0.056 

2.340 

16. 

Ultimate  strength 

57  ksi 

114.22 

17. 

a:  Shape  parameter 

4.0 

35.59 

18. 

fj : Scale  parameter 

30,000  hrs 

15.37 

19. 

U-j_:  Probability  of  inspecting  cracked 
details 

1.0 

129.13 

20. 

a^:  Maximum  undetectable  crack  size 

0.02  in 

0.457 

21. 

a^:  Minimum  crack  size  definitely  detectable 

0.3  in 

2.002 

22. 

m:  Parameter  appearing  in  expression  for 
U2(a) 

0.125 

8.242 

23  . 

(AI'k)l/k:  Threshold  value 

1.5  ksi  /in 

532.11 

24  . 

X 1 : Slope  of  crack  propagation  law 

10.0 

152.73 

25. 

X : Slope  of  crack  propagation  law 

2.5 

54.56 
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TABLE  6. 


REPRESENTATIVE  VALUES  OF  SHAPE  PARAMETER  a 


(a)  Between  10  and  20  percentile 


Material 

Fighter 

Transport 

Aluminum 

4.5 

3.5 

Titanium 

3.0 

2.5 

Steel  (100-200  ksi) 

3.5 

3.0 

Steel  (200-300  ksi) 

2.5 

2.0 

(b)  In 

the  median  range 

Material 

Fighter 

Transport 

Aluminum 

3.0 

6.0 

Titanium 

6.5 

4.0 

Steel  (100-200  ksi) 

7.0 

5.0 

Steel  (200-300  ksi) 

5.0 

3.5 
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TABLE  7.  SCATTER  FACTORS  S FOR 
RELIABILITY  LEVEL  R = 0.5 

(a)  n = 1 


2 

3 

4 

5 

3 

1.7 

1.4 

1.3 

1.3 

25 

5.0 

2.9 

2.2 

1.9 

100 

10.0 

4.6 

3.2 

2.5 

250 

15.8 

6.3 

4.0 

3.0 

1000 

31.6 

10.0 

5.6 

4.0 

(b)  n = 3 

m\a 

2 

3 

4 

5 

3 

1.9 

1.5 

1.4 

1.3 

25 

5.6 

3.2 

2.3 

2.0 

100 

11.3 

5.0 

3.4 

2.6 

250 

18.0 

6.9 

4.2 

3.2 

1000 

36.8 

10.9 

5.9 

4.2 
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TABLE  8.  SCATTER  FACTORS  S FOR 


RELIABILITY  LEVEL  R = 0.5 


(a)  n = 1 


m\a 

3 

2 

5.2 

3 

3.0 

4 

2.3 

5 

1.9 

25 

15.0 

6.1 

3.9 

3.0 

100 

30.0 

9.7 

5.5 

3.9 

250 

47.4 

13.1 

6.9 

4.7 

1000 

94.9 

20.8 

9.7 

6.2 

(b)  n = 3 

2 

3 

4 

5 

3 

5.3 

3.0 

2.3 

1.9 

25 

15.3 

6.2 

3.9 

3.0 

100 

30.6 

9.8 

5.6 

3.9 

250 

48.3 

13.2 

7.0 

4.8 

1000 

97.1 

21.0 

9.8 

6.3 
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TABLE  9.  SCATTER  FACTORS  S FOR 
RELIABILITY  LEVEL  R = 0.99 


(a)  n = 1 


m\a 

2 

3 

4 

5 

3 

17.2 

6.7 

4.2 

3.1 

25 

49.8 

13.5 

7.1 

4.8 

100 

99.5 

21.5 

10.0 

6.3 

250 

157 

29.1 

12.5 

7.6 

1000 

314 

46.3 

17.7 

10.0 

(b)  n = 3 

m\a 

2 

3 

4 

5 

3 

17.0 

6.7 

4.2 

3.1 

25 

50.0 

13.5 

7.1 

4.8 

100 

99.0 

21.5 

10.0 

6.3 

250 

160 

29.1 

12.5 

7.6 

1000 

310 

46.3 

17.7 

10.0 
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